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Preface 



This short lecture is consist of three parts. 

In the first part, the Liischer's formula, which relates the scattering phase shift 
to the two particle energy in the finite volume, explained. A comprehensive but less 
rigorous derivation for the formula has been attempted in this lecture for the tttt 
system as an example with the emphasis on the Bethe-Salpeter(BS) wave function. It 
is important to note that the BS wave function at large separation behaves as the free 
scattering wave with the phase shift which is determined by the unitarity of the S— 
matrix in QCD. The Liischer's formula can be obtained from this asymptotic behavior 
of the BS wave function. 

In the second part, the BS wave function is considered at non- asymptotic region 
where the interaction between two particles exists, in order to define the potential in 
quantum field theories. This method is applied to the two-nucleon system in order to 
extract the NN potential from lattice QCD. 

In the last part, the origin of the strong repulsion at short distance in the NN po- 
tential, called the repulsive core, is theoretically investigated by the Operator Product 
Expansion and the renormalization group. 



Acknowledgements 



I would like to thank the organizer of this Les Houches school for giving me an op- 
portunity for this lecture at the school and all attendees of the school for stimulating 
discussions. I also thank members of HAL QCD(Hadron to Atomic nuclei from Lattice 
QCD) Collaboration, T. Doi, T. Hatsuda, Y. Ikeda T. Inoue, N. Ishii, K. Murano, K. 
Nemura and K. Sasaki, for providing me the latest data and useful discussions, and 
my collaborators and friends, J. Balog, N. Ishizuka, W. Weise and P. Weisz, for valu- 
able discussions and comments. This work is supported in part by Grant-in- Aid of the 
Japanese Ministry of Education, Sciences and Technology, Sports and Culture (Nos. 
20340047) and by Grant-in- Aid for Scientific Research on Innovative Areas(No. 2004: 
20105001, 20105003). 



Contents 



1 Intorduction: Nuclear Forces 

2 Phase Shift from Lattice QCD: Liischer's formula in the finite 
volume 

2.1 Bcthe-Salpeter wave function and phase shift in quantum field 
theories 

2.2 Liischer's formula for the phase shift in the finite volume 

2.3 Some references for the tttt phase shift from lattice QCD 

3 Nuclear Potential from Lattice QCD 

3.1 Strategy to extract potentials in quantum field theories 

3.2 Extraction of the BS wave function on the lattice 

3.3 Tensor potential 

3.4 Results in lattice QCD 

4 Repulsive core and operator product expansion in QCD 

4.1 Basic idea 

4.2 Rcnormalization group analysis and operator product expansion 

4.3 OPE and Anomalous dimensions for two nuclcons 

4.4 Short distance behavior of the potentials and the repulsive core 

5 Concluding Remarks 



References 



1 

Intorduction: Nuclear Forces 



In 1935 Yukawa introduced virtual particles, pions, to explain the nuclear force( Yukawa, 1935 ) 
which bounds protons and neutrons inside nuclei. Since then, enormous efforts have 
been devoted to understand the nucleon-nucleon (NN) interaction at low energies both 
from theoretical and experimental points of view. To describe the elastic nucleon- 
nucleon (NN) scattering at low-cnergics below the pion production threshold to- 
gether with the deuteron properties, the notion of the NN potential turns out to 
be very usefulCTakctani et ai, 1967 Hoshizaki et al, 1968 Brown and Jackson, 1976 



[Machleidt, 1989. ; >lachleidt and Slaus, 2001D : it can be determined phenomenologi- 
cally to reproduce the scattering phase shifts and bound state properties. Once the 
potential is determined, it can be used to study systems with more than 2 nucleons 
by using various many-body techniques. 

Phenomenological NN potentials which can fit the TV TV data precisely (e.g. more 
than 2000 data points with x^/dof ~ 1) at Tiab < 300 MeV are called the high- 
precision TV TV potentials. Those in the coordinate space, some of which are shown in 
Fig ll.li are known to reflect some characteristic features of the TV TV interaction at dif- 



ferent length scales ( Taketani et ai, 1967 Hoshizaki et ai. 
[Machleidt, 19891 [Machleidt and Slaus, 2001D : 



1968 Brown and Jackson, 1976 



(i) 



(ii) 



(iii) 



The long range part of the nuclear force (the relative distance r > 2 fm) is domi- 
nated by the one-pion exchange introduced by Yukawa ( Yukawa, 1935[ ). Because 
of the pion's Nambu-Goldstone character, it couples to the spin-isospin density of 
the nucleon and hence leads to a strong spin-isospin dependent force, namely the 
tensor force. 

The medium range part (1 fm < r < 2 fm) receives significant contributions from 
the exchange of two-pions (tttt) and heavy mesons (p, lo, and a). In particular, the 
spin-isospin independent attraction of about 50 - 100 MeV in this region plays 
an essential role for the binding of atomic nuclei. 

The short range part (r < 1 fm) is best described by a strong repulsive core 
as originally introduced by Jastrow ( [Jastrow, 1951[ ). Such a short range repul- 
sion is important for the stability of atomic nuclei against collapse, for deter- 
mining the maximum mass of neutron stars, and for igniting the Type II super- 



20001 



nova explosions ( Tamagaki et al, 1993[ [Heiselberg and V. Pandharipande, 2000 
Lattimer and Prakash 



A repulsive core surrounded by an attractive well is in fact a common feature of the 
"effective" potential between composite particles. The Lenard- Jones potential between 
neutral atoms or molecules is a well-known example in atomic physics. The potential 
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Fig. 1.1 Three examples of the modern A'^A'' potential in ^So (spin singlet and 
S-wave) channel: CD-Bonn (Machleidt, 200 Reid93( [Stoks et al, 1994[ ) and Argonne 
fi8( |Wiringa et al, 1995| ). 

between ^He nuclei is a typical example in nuclear physics. The origin of the repulsive 
cores in these examples is known to be the Pauli exclusion among electrons or among 
nucleons. The same idea, however, is not directly applicable to the NN potential, 
because the quark has not only spin and flavor but also color which allows six quarks 
to occupy the same state without violating the Pauli principle. Therefore, to account 
for the repulsive core of the NN force, various ideas have been proposed as sum- 
marized in Ref. ( |Myhrer and Wroldsen, 1988[[Oka et al, 2000l [Fujiwara et al, 2007D : 
exchange of the neutral uj meson (|Nambu, 1957[), e xchange of non-linear pion field 
( Jackson and Pasquier, 1985 Yabu and Ando, 1985 ), and a combination of the Pauh 
principle with the one-gluon-exchange between quarks (Otsuki et al, 1965!'Machida and M. Namiki, 1965 
Neudachin et al, 1977||Liberman, 1977||DeTar, 1979||Oka and Yazaki, 1980^ Oka and Yazaki, 1981a| 



Oka and Yazaki, 1981 6[ Toki, 1980{ Faessler et al, 1982 ). Despite all these efforts, con- 
vincing account of the nuclear force has not yet been obtained. 

In this situation, it is highly desirable to study the TV TV interactions from the first 
principle lattice QCD simulations. A theoretical framework suitable for such purpose 
was first proposed by Liischcr ( [Liischer, 1991 ): For two hadrons in a finite box with 
the size L x L x L in periodic boundary conditions, an exact relation between the 
energy spectra in the box and the elastic scattering phase shift at these energies 
was derived: If the range of the hadron interaction R is sufficiently smaller than the 
size of the box R < L/2, the behavior of the two-particle Bethe-Salpeter (BS) wave 
function -(/"(r) in the interval R < |r| < L/2 under the periodic boundary conditions 
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has sufficient information to rciatc the phase shift and the two-particle spectrum. 
Liischer's method bypasses the difficulty to treat the real-time scattering process on 
the Euclidean lattice. Furthermore, it utilizes the finiteness of the lattice box effectively 
to extract the information of the on-shell scattering matrix and the phase shift. 

Recently, a closely related but an alternative approach to the NN interactions from 
lattice QCD has been proposed jlshii et al, 2007l|Aoki et al, 2008|[Aoki et al, 2010&D . 
The starting point is the same BS wave function ■0(r) as discussed in Ref. ( Liischer, 1991[ ). 
Instead of looking at the wave function outside the range of the interaction, the au- 
thors consider the internal region |r| < i? and define the energy- independent non-local 
potential U{r, r') from ip{r) so that it obeys the Schrodinger type equation in a finite 
box. Since U{r,r') for strong interaction is localized in its spatial coordinates due 
to confinement of quarks and gluons, the potential receives finite volume effect only 
weakly in a large box. Therefore, once U is determined and is appropriately extrapo- 
lated to L ^ OD, one may simply use the Schrodinger equation in the infinite space to 
calculate the scattering phase shifts and bound state spectra to compare with experi- 
mental data. Further advantage of utilizing the potential is that it would be a smooth 
function of the quark masses so that it is relatively easy to handle on the lattice. This 
is in sharp contrast to the scattering length which shows a singular behavior around 
the quark mass corresponding to the formation of the NN bound state. 

In this lecture, we first introduce the Liischer's method for the scattering phase 
shift in SecEl Since the method is not only well established but also well explained in 
Ref. ( [Liischer, 1991] ), we mainly consider properties of the BS wave function, in terms of 
which the scattering phase shift can be related to the energy shift of the 2 particles state 
in the finite box. These properties are also used to define the NN potential in SeclSl 
where new method in Ref.( |Ishii et al, 2007| |Aoki et al, 2008[ |Aoki et al, 20106[ ) is 



explained in detail. We finally consider a very recent attempt to understand the origin 
of the repulsive core in the NN potential in SecH) Using the operator product expan- 
sion and the renormalization group analysis in QCD, the potential derived from the 
BS wave function in Secl3]is shown to have the repulsive core, whose functional form 
is also theoretically predicted(?). Brief concluding remarks are given in SecH] 
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Phase Shift from Lattice QCD: 
Luscher's formula in the finite 
volume 



2.0.1 Preparation: Scattering phase shift in quantum mechanics 



In this subsection, as a preparation for the latter sections, we give some basics of the 
scattering theory in quantum mechanics. 

Let us consider the 3-dimcnsional Schrodinger equation, given by 



[Ho + V{v)]^{v) = E^{v) 



where 



2m' 



(2.1) 



(2.2) 



Hereafter we consider only the spherically symmetric potential that ^(r) = V{r) with 
r = |r|. In this case it is convenient to use the polar coordinate such that 

1 



df dr ' 



1 d 



sin 6 I 



sin 9 



d_ 

do 



sin^ 6 ' 



(2.3) 



Using the separation of the variables, we consider the following form of the solution 
V3(r): 



(2.4) 



where the spherical harmonic function satisfies 



and is normalized as 

r2n 



sin Ode Yim{0A)Yi'. 



Sii/S, 



(2.5) 
(2.6) 

(2.7) 



Note that in this lectrue X means a complex conjugate of X while X'' is a hermitan 
conjugate of X. Explicitly it is given by 



'2/ + ! il-\m\y. 
An ' {I + \m\)l 



Phn{cos9)e 



im4> 



(2.8) 
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(-1)'" m > 
1 m < 



where Pim{z) is the Legendre bi-function of degree /, defined by 



Pi, 



(2.9) 

(2.10) 
(2.11) 



For small / and m, for example, we have 
1 



5^20 ( 

Y2M 



47r 
5" 1 



■0086*, yi±i = =Fa/— sin^e^^"^ 
V Stt 



15 



— -(3cos^6'- 1), y2±i(6i,0) = sin 61 cos 6*6 



±i4 



47r2 



Stt 



15 
32^ 



sin^±'2,^. 



(2.12) 



From eqs. <\2A\ and (|2.5p . the 3-dimensional Schrodinger equation (|2.ip is reduced 
to the 1 dimensional equation for Ri as 



(rRdr)) + (2m (i? - F(r)) - 



i?i(^) -0. 



Usually we further assume the following properties for the potential V{r) 

lim rV(r) = 0, 



lim r"V{r) = for V?! G Z. 



The first condition means 



1 



for small r, which leads to 

Ri{r) = 0{r') 
for small r. The second condition means 

V{r) =0 for r > i?, 
for sufficiently large R, so that eq. (|2.13p becomes 



^ 2 , r ?(? + l) 
i?; (y) + -R{y) + \l- ±^ 

y I y 



Ri{y)^o 



(2.13) 



(2.14) 
(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 
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for r > R, where y = kr and fc^ ~ 2mE. The solution to this equation is obtained as 

Riiy) ^ Aijiiy) + Bmiiy) (2.20) 
where spherical Besscl functions are given by 



. N . Ml dV famx\ (2Z + 1)!!' 
Jiix)^{-xy —] —]^{ , (2.21) 



X dx 



sm(x — Itt /2) 

^ X oo 

X 

x^O 



ni{x) = {-xy[-—] . (2.22) 



X dx J \ X 



cos{x — Itt/2) 



X 

Therefore, in r — > oo limit, the above solution becomes 

, , sin (fcr — ^tt) cos (fcr — ^tt) sin (fcr — ^tt + (5;(fc)) , 

Ri{r)^Ai ^ ^—L + Bi ^ ^ f ^ 2.23 

kr kr kr 

where the scattering phase shift 5i{k) is given by 



tan 8i{k)^^ , Ci = ^Aj+ Bf (2.24) 



and Ai — Ci cos6i{k), Bi = Cism5i{k). If 5i{k) > the interaction is attractive for 
this k, while it is repulsive if Si{k) < 0. 

For the particle-particle scattering in quantum mechanics, we consider the La- 
grangean, 

L = ^mirj + ^m.Jl - V{\ri - r2|). (2.25) 

By introducing the relative coordinate r = ri — r2 and the center of gravity R = 
(miTi + 7112^2) / {mi + "^2)- the above Lagrangean becomes 

L = ^MR? + i^f2 - V{r), (2.26) 

where M = mi + m2 is the total mass and — mim2 / {mi +7712) is the reduced mass. 
The corresponding Hamiltonian is given by 

H^Hg + i/rei, Hg = i^.ci = ^p' + V{r), (2.27) 

where P = MR and p = /ir. While Hg is a Hamiltonian for a free particle, H^ei corre- 
sponds to the Hamiltonian for a particle under the potential V{r), whose Schrodinger 
equation identical to eq. ()2.ip . 
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2.1 Bethe-Salpeter wave function and phase shift in quantum 
field theories 

In this section, wc construct a "scattering wave" in quantum field theories whose 
asymptotic behavior is identical to the one in the scattering wave in quantum me- 
chanics given in ea. (|2.23p . Moreover we show that the phase shift corresponds to the 
phase of the S'-matrix required by the unitarity. For notational simplicity we consider 
the TTTT scattering in QCD here. 

The unitarity of the S'-matrix = SS^ = 1 with S = 1 + iT leads to 

(/|T|z) - (/|rt|z) ^ ^J2{f\T^\n){n\T\^), (2.28) 

n 

where \n) are asymptotic states. In the case of tttt scattering in the center of mass 
system that ka + ki, —i' kc + kd where ka = (-Efc, k), k^ = [Ek, — k) and k^ = (i?p, p), 
kd = {Ep, — p) with Ek = \/k^ -|- and Ep = ^p^ -I- mj, we explicitly write 

(fce,fcd|T|fc„fc6) = (27r)4,5W(fc, + fcb-fc,-fcd)r(p,q) (2.29) 

We consider the elastic scattering, where the total energy is below the 47r produc- 
tion such that 2\JkP- 4- < 4m7r, equivalently, fc^ < 3m^ with = |k| = |p|. In this 
case, due to the energy-momentum conservation, the sum over intermediate state n in 
Eq. (|2.28p is restricted to the tttt states as 

Inserting this into Eq. (|2.28p . we have 

r(p,k)-rt(p,k) = ^ / dr!qTt(p,q)r(q,k), (2.31) 

where |q| = A: and Slq is the solid angle of the vector q. Using the partial wave 
decomposition that 

oo / 

T(p,k) = 4^^ Ti{k)Yir^{n^)Yim{Sl^) (2.32) 

;=o ■ni=-l 

and the orthogonal property of the spherical harmonics function y;„j(0, (j)) ~ i^,ri(^^q) 
that 

J dflqYlmi^q) Yl'm'i^q) = Sll'S„im' , (2.33) 

the unitarity (|2.31|) becomes 

Ti{k) -Tdk) = i-^W^ Ti{k). (2.34) 
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A solution to this unitarity condition is easily obtained as 

T,(fc) = l^e'*'Wsin<5,(fc), (2.35) 
k 

where Si{k) is an arbitrary real function of fc, and can be interpreted as the scattering 
phase shift, as seen later. 

Wc now introduce the Bethe-Salpeter (BS) wave function for tttt system, defined 

by 

ifir) ^ (0|r{^a(x + r, ta)TTb{^, tb)}\ka,a, kb, b; in) (2.36) 

where \ka,a,kb,b;m) is a tttt asymptotic in-state in the center of mass system such 
that ka = {Ek,^) and kb = (-E^,— k) with flavors a and b. The pion interpolating 
operator is given by 

7ra(x,xo) = q{x)i-f5Taq{x), q{x) = (^^[j^j^ (2-37) 

with the Pauli matrix Ta . For simplicity we take ta ~ tb + e with e > and the e — > 
limit. Wc then simply write 

ifiir) = {0\TTaix + r,t)TTb{^,t)\ka,a,kb,b;m). (2.38) 

The name, the Bethc-Salpeter wave function, comes from the fact that this quantity 
satisfies the Bethe-Salpeter equation! Bethe and Salpeter, 1951 ). Unlike field equations 
such as the Dyson-Schwinger equations, the BS equation is derived from the diagram- 
matic consideration. We here consider the BS wave function from the different point 
of view. 

By inserting the complete set of the out-states such that 

/■ d^p , , , , V- |^:out)(X;out| 



we have 

ip{v) = (^'='^"'"=(r) + <^i"°ia««<=(r) (2.40) 

where 

^eiastic(^) = / "^'j e-P-(x+r)-^pot(p^ ^. out|^b(x, ^)|k, a, -k, b; in) 

J (27r)-^2po 

(2.41) 

^inelastic ( ^ ^(0|7r,(x + r, t)\X; out)-^{X; out|7rb(x, t)\k, a, -k, b; in) (2.42) 
and |p,a, ;out) is an one-pion out-state with momentum p, which satisfies 



(OKa(x,xo)|p,6;out) = SabVZ^e-'P-^ , po - ^p^+ml- (2.43) 
On the other hand, X represents general states other than one-pion states. 
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Using the reduction formula that 

aout(p)r(0) - T(0)ain(p) = (V + ml)T{Tr{p)0) (2.44) 

T{0)al{p) - al,{p)T{0) = (V + ml)T{OnHp)) (2.45) 
where O is an arbitrary field operator and 

/ ipx 

d*x-^n{x), (2.46) 

we obtain 

(p,a;out|7r6(x,t)|k,a,-k,6;in) = ^/z^{2Trf2koS^^\p - -k)e-''''' 



Z^— ^ -T{p,q,ka,h) (2.47) 

where off-shell T-matrix T is defined by 

f{p, q, ka, kb) ^ {-p^ + ml){-q'^ + mJ)G(p, q, ka, kb){~kl + ml){~kb + m^) 

(2.48) 

G{p, q, fc,.fcb)i(27r)4^(4)(p + q^ka^kb)^ (0|T{7r,(p)7rb(g)4(A:a)4(fcb)}|0). 

(2.49) 

Here p = {p,po), ka = (k, fco) and kf, = (— k, fcg) are on-shell 4-momenta, while 
q = (— p, 2fco — Po) is generally off-shell. Using this expression we obtain 

TTT^ 2 ^—T{p,q, ka,h) 

(27r)'^2po mi~ q^ ~ le 

(2.50) 

Similarly we have 



inelastic/„\ -i2fcot \J Z-j^Zx e'^"'^'' - , ; ; \ /o r;i\ 

<^ r) = e Z^^TT^ 5 2 -Tx[px,q,ka,kb) (2.5f) 

^ 2i/x — q^ — le 

where q — (— pjc, 2fco — (px)o)- For simplicity we hereafter set t = 0. We rescale (^'=1'^^'"= 
as 

(^(r) = Z^^^''^^*''=(r) + <pi"<=ia''ti'=(r) (2.52) 

where 

^eia.tic(^) = e^"-- + / ^-^-±--i/(p,k)e-- (2.53) 



with 



and k — |k| and p = |p| 



f/(p,k) = ^±^r(p,q), (2.54) 
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We now investigate the large r = |r| behavior of the BS wave function below 
the An inelastic threshold. We first consider (^<=iastic^ Using the following partial wave 
decomposition 



if(p,k) = 'i7rJ2Hiip,m„,{np)Yi^ini,) (2.55) 

^olastic(j.) ^ 4^^i'^clastic^^^ fc)F/„ (f^OW^O (2-56) 

e»P- = 47rY,i' 3l{pr)Yln^{n,)Ylm{np) , (2.57) 



we have 



^f-«^(r, k) = Mkr) + I l^^-^—H,{p, kMpr). (2.58) 
We assume that the interaction vanishes for large r: 

- (V2 + fc2)^eiastic(^) ^ I ^H{p, k)e^P- — , (2.59) 



which, in terms of the partial wave, gives 

p^dp 



27r2 



Hi{p,k)ji{pr) ^ 0. (2.60) 



We evaluate the second term of cq. p.58p . Using the explicit form of the spherical 
Bessel function ji{x) in cq. p.2ip we have 



27r2 p'^ ^ k^ — ie ' dr J Jq An'^ipr p"^ — k^ — ie 

1 dV f°° p^-^dp 



^ \r dr ) J _^ An'^ipr p'^ — k'^ ~ ie 

(2.61) 

Here we used the property that Hi{—p, k) — {—lY Hi{p, k), which is shown as 
H{p, k) = 4TrY,Hiip, fc)yi™(r!p)y/„(r!k) = 47r^ i//(-p, k)Yi„,in-p)YhJn^ 

l,m Urn 

= 47r^(-i)'i/,(-p, fc)yi™(rip)y,„(r!k) (2.62) 

where fi^p means 6 ^ n — 9 and — > tt + 0. To proceed, wc assume Hi{p, k) — 0{p^) 
for small p, so that no contribution around p = Q appears in the above integral. 

We know that the half off-shell T-matrix Hi (p, k) do not have any poles and cuts in 
the real axis, since k"^ is smaller than inelastic threshold. For simplicity we assume that 
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Hi(j), k) has only poles in the complex p plain(this argument may also be generalized 
for cuts): 



E 



E 



n>OJmA;„>0 n<0,ImA;„<0 



(2.63) 



where Hi(j), k) is analytic in p. Note that fc„ and Z„ implicitly depend on k. This form 
of the assumption satisfies the condition eq. p.GOp . which becomes 



1 d 
r dr 



E 



n>0,lra/c„>0 



" " ^»(Rcfcn)r^-(Imfc„)r ^ q 

27rfc„r ' 



(2.64) 



for /co^" 2> 1, where fco = min,i>oImfc,i > 0. Using the assumption (j2.63p . we can 
evaluate ea. p.6ip as 



47r 



Hi{k,k){~krf 



1 d 



kr dkr J kr 



^ ^ikr 



+ i-ry 



1 d 



1 d 



r dr J ^ 27rfc„r(fc2 _ fc2) 



—H^ik, k){ni{kr) + zjiikr)} + (-r)' (-j;;) 



r dr 



n>0 



l-Kknrikl - P) 



(2.65) 



Since Imfc„ > 0, the sum over n > vanishes exponentially for large r which satisfies 
fco?" ^ 1- Similarly we can show that (^mciastic j-j.^ vanishes exponentially for large r as 
long as k^ is smaller than inelastic threshold. 

Therefore, for large r (kor ^ 1), we finally obtain 



Ji{kr) + —Hi{k,k){ni{kr) +iji{kr)} 
An 



JSi(k) 

kr 



■sin(fcr + <5/(fc)-/7r/2) 



(2.66) 



where the folowing form of the on-shell T- matrix determined from unitarity in eq. (|2.35p 



4-7r 

H(k,k) = — e^-^'W sinSiik), 
k 

and the asymptotic behaviour of ji{x) and ni{x) that 

sin(x — /7r/2) cos(x — /7r/2) 
ji{x) ~ , ni{x) ~ 



(2.67) 



(2.68) 



are used. 

Note that the derivation for the large r behavior of the BS wave function in this 
section is similar but a little different from that in (Lin et ai, 2002 Aoki et ai, 200561 
for TTTT and in ( [Ishizuka, 2009| |Aoki et ai, 2010"6| for NN, though the final results are 
same. 
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2.2 Luscher's formula for the phase shift in the finite volume 

We now consider the finite volume fLiischer, 1991 1. We assume that no interaction 
(except exponentially small contributions) exists at r > i?, where R is sufficiently large. 
Therefore, if the box size L is larger than 2R, there exist a region that R < r < L/2 
where 

(V^ + A:^)(^L(r;fc) = (2.69) 

is satisfied for the BS wave- function V'L(r; k), which is given by 

</7L(r;/c) = (0|^a(x + r,0)7rb(y,0)|fca,a,fc6,6)L, (2.70) 

where the subscript L indicates that the state is constructed in the finite box. This 
wave function is expanded in terms of the BS wave function in the infinite volume, 
introduced in the previous sections as 

M-^;k) = 'inJ2CiMMr,k)Yi,ninr) (2.71) 

l,m 

where the coefficient Cim{k) is introduced to satisfy the periodic boundary condition 
such that <^L(r + nL; k) = (/^^(r; fc) for n — {nx, Uy, Uz) G Z"^. Note that 

ipiir,k) ^ ni{kr)e'^'^''^ sinSiik) + ji{kr)eJ'^''^''^ cos 5i{k) (2.72) 

for r > R. 

On the other hand, we can construct the solution of the Helmholtz equation 

(V2 + fc2)^i(r;fc) = (2.73) 
for r 7^ with the periodic boundary condition as 

Mr; k)=Y, vim{k)Gim{r, k) {2.1 A) 

where 

Gi™(r, k) = %/4^Y,„(V)G(r, k) (2.75) 

G(r''=)==T¥E^CT' r=|p|p = n^,nGz4, (2.76) 
per v \. ) 

Yz„,(p) =p'yz„,(r!p), p=|p|. (2.77) 

It is easy to see the above (pL(r;fc) satisfies both Helmholtz equation and periodic 
boundary condition for arbitrary vim{kys as 

(V^ + fc2)^i(r; k)^Y. Vhn{k)^/A^Yi^{\/){y^ + e)G{v, k) 
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gi,n{k) = V^— J2 p2 _fc2 ^''"("p)- (2-81) 



= ^i'/™(fc)%/4^Yz„,(V),5(3)(r) = (2.78) 

Z,m 

for r 7^ and 

(/3L(r + nL;fc) = ^ Wi,„(fc)V4^Y/,„(V)G(r + nL, fc) 

l,m 

= ^t./™(A:)V4^Yz„,(V)G(r,fc) = (^L(r;fc) (2.79) 

The coefficient vim can be determined by comparing cq. ()2.74p with eg. p. 71 p . We 
first rewrite 

G(r, k) = -^noikr) + V VI^Yira{n,)gimik)ji{kr) (2.80) 
47r 

where 

1^ {ip/kY 

per 

This can be easily seen as follows. Since 

(V2 + fc2)A^„(fcr) =,5(r), (2.82) 
47r 



G(r,fc) - — no(fcr) (2.83) 
47r 

satisfies the Helmholtz equation for all r and is smooth at r ^ 0, so that it can be 
expanded by ji as 

G(r,fc) - ^no(fcr) = ^ V4^<7ir„(fc)j,(fcr)l^,„(r!r). (2.84) 

Lrn 

Using 

e'P-- = 4^ ^ z'iz (pr)l^„ (Or)rz™(r!p) (2.85) 

lrn 

in cq. (|2.76p . and considering the r ^ limit in the both side of eq. (|2.84p . we obtain 
eq.dlHI). (Note tha t ji{x) - x'-/{2l + 1)!! as a; ^ 0.) 
We next observe ( pjiischer, 1991D that 

Gim{r,k) = VI^YimiV)Gir,k) = t^lh 

47r 
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Table 2.1 Non-zero independent elements of 





a (Afoo) 


b ( A'/40 ) 


c (A/eo) 


Moo, 00 


1 










1 








M20,20 


1 


18/7 





M21,21 


1 


-12/7 





M22,22 


1 


3/7 





M22,2-2 





15/7 





Mso^io 





-4V^/7 





M3i,ii 





3\/T4/7 





M33,i_i 





v^/7 





Mao^ao 


1 


18/11 


100/33 


Mai,ai 


1 


3/11 


-25/11 


M32,32 


1 


-21/11 


10/11 


Ma2,3_2 





15/11 


-70/11 


Maa^aa 


1 


9/11 


-5/33 


Maa,3-i 





3^15/11 


35VT5/33 



m,l'ra' 

I' ,m' 

where non-zero elements of Mim^i'm' arc given by the linear eombination of 

1 47r 



Mir, 



-gim{k). 



i\2l + l) k 
The following properties generally hold: 

Non-zero elements at /, /' < 3 are expressed as 



(2.86) 



(2.87) 



(2.88) 



(2.89) 



with a, 6, c given in table 12. ll See Ref.( |Luscher, 1991 1 for more details. 

We now consider the cubic group 0(3, Z), which has 24 elements and is generated 
by following elements in the special cubic group SO (3, Z) 



1 

i?x = I -1 I , ii-j, = 
1 




(2.90) 



and the parity transformation Pr = — r. There are five irreducible representations 
in S0(3,Z), denoted by Ax, A2, E, Ti and T2, whose dimensions are 1,1,2,3, and 
3, respectively. The irreducible representations of 0(3, Z) are constructed these five 
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Table 2.2 Decomposition of the angular momentum into irreducible representa- 
tions of the cubic group 



I 


rep. 


basis polynomials 


independent elements 





At 


1 




1 






i = 1,2,3 


2 


E+ 


2 2 


(i,j) = (l,2),(2,3) 


2 






i 7^ 3 


3 








3 






I = 1,2,3 


3 


T2 




(i,j,fc) = (l,2,3),(2,3,l),(3,l,2) 


4 


At 


5(rt + r|+r|)-3r4 




4 


E+ 


7(^4 4) _ 6^2(^2 _ ^2) 


(z,j) = (l,2),(2,3) 


4 




r.ir^ — rjrf 




4 




7{rir^ + rjrf) — dr^rirj 


i 7^ 3 



irreducible representations of S0(3,Z) and the parity eigenvalue ±1. It is noted that 
irreducible representations of the rotational group 0(3,R) are decomposed in terms 
of these irreducible representations. For example, 

= A+, 1 = Ti", 2 = E+ (S T+, 

3 ® Tf ® T^, 4 = A+®E+®T+® T+ , (2.91) 

where a number denotes an eigenvalue of the angular momentum I. The corresponding 
basis polynomials for each cubic representation are given in tabl d2.2l 

We now compare two expressions of ipL(j',k) in some irreducible representations 
of the cubic group. If we project the BS wave function to representation, which 
contains the / = partial wave as well as Z > 4 contribution. Neglecting I > 4 
contribution, ea. (|2.71l) becomes 

ip'^'{r;k) = V4^Coo(fe)e"^''^''' [no(fcr) sin (5o(fc) + jo{kr) cos5o{k)] (2.92) 
for r > R. In order to match this expression, eq. ((2J4)) must be 



(/3^' {v\k) = V47rwoo(fc)G'oo(r, fc) 



47rW( 



'00 



no(fcr) + ^m.imXim{^v)3iikr) 



Im 



(2.93) 



since G;„i(r, k) with / 7^ 0, which contains ni{kr), can not appear in this equation. By 
comparing the two, we have 

Coo(fc)e*'°(")sin5o(fc) =«oo (2.94) 
Coo(fc)e''^«(*'-) cos5o(fc) = «ooMoo,oo = «ooMoo, (2.95) 

which leads to the famous Liischer's formula ( [Liischer, 1991[ ), 
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cot{So{k)) = Moo = y.goo(fc) = E (2.96) 

per ^ 

Note that un-matched components proportional to Moo, 4m ("^ = 0,±4) give / = 4 
contributions. 

Let us briefly explain how to use this formula. We first calculate the energy E2{L) 
of two pions in the center of mass frame on the finite box with the periodic bound- 
ary condition, where L is assumed to be larger than 2R. We then determine k, the 
magnitude of the relative momentum of the two pions, from the equation that 

E2{L) = 2v/fc2 + m2, (2.97) 

where ttitt is the pion mass in the infinite volume limit. We finally determine Sa{k), 
by solving eq. ()2.96|) with k from eq. ()2.97p and the spatial size L. It should be noted 
that the momentum for one pion is quantized as p = 2T:n/L on the finite box with 
the periodic boundary condition. If an interaction between two pions were absent, we 
would have E2{L) = 2-\/p2 + m^, so that k = |p|. The presence of the interaction 
makes fc is a little different from |p| on the finite box. The above Liischer's formula 
relates the difference to the scattering phase shift. By applying the formula for two 
pions with different p and L, we can determine the scattering phase shift 5o{k) at 
several values of k. 

Similarly as above, we consider the Tj~ representation, which contains the I = 1 
partial wave as well as Z > 3 contributions. If latter ones are neglected, we again obtain 

cot(<5i(fc)) = A/oo--^5oo(fc). (2.98) 

On the other hand, if the contribution from I = 3 partial wave can not be neglected, 
in addition to the / = 1 component, we have 

/i}{v-k) = 47r^Cim(fc)yim(17r.)e^^i('=) [ni(fcr) sin5i(fc) + ji(A:r) cos^i(fc)] 

m 

+ 47r^C3m(fc)r3m(^^r)e'*^''=^ [rig (fcr ) sin .^g (A:) + J3 (^0 COS (Jg (fc)] , 
m 

(2.99) 

which should be compared with 

/L^{Y-k) = Y,^,^^n{k)Gi,n{y,k) + Y,v^m{k)G:i^^{v,k). (2.100) 
From the matching condition we obtain, after a little algebra, 

d^tf^'^M "''^'^'^'m ^"'".Aml-O' (2-101) 
M3o,io, M3o,3o - cotd3(fc)y 

for m = 0, and 
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/Mii,ii -cot(5i(fc), M31.11, M33,i_i 
det M3i,n, M3i,3i -cot,53(fc), M33,3_i | = 0, (2.102) 

\ M33,i_i, M33^3_i, M33^33 - COt(53(fc) 

for m = 1, —3 or m = —1, 3, while we have 

det f - cot^3(fc), M3,3 . \ ^ (2^03) 

\^ M32,3-2, M32,32 - COt(53(fc)y 

for m — 2,-2. From the last equation, we can determine (53(fc) at some k and L. 
Putting this S3{k) into the first or the second equation, we can also extract Si{k). 

2.3 Some references for the tttt phase shift from lattice QCD 

The BS wave function for the tttt system in the isospin 1 = 2 channel has been 



investigated in quenched QCD at fc ~ 0(Aoki et ai, 2005a) to extract the scattering 



length flo through the Liischer's formula in the center of mass system. The scattering 
length is related to the scattering phase shift So{k) as 

-4m=- + ^0T + O(^')' (2.104) 

tanoo(fcj flQ z 

where tq is called the effective range. The calculation of the BS wave function for 
the TTTT system in the 1 = 2 channel has been extended to the case of the non-zero 
momentum in quenched QCD( [Sasaki and Ishizuka, 2008[ ) and the / = 2 tttt scat- 
tering phase shift can be extracted using the Liischer's formula in the laboratory 
system( Rummukainen and Steven Gottlieb, 1995[ ) . 



Besides quenched calculations, there are only a few calculations for the tttt scatter- 
ing phase shift. 

The / = 2 TTTT scattering length and phase shift have been calculated through the 
Liischer's formula in 2-flavor lattice QCD with the 0{a) improved Wilson fermion in 
both center of mass and laboratory systems ( jYamazaki et al., 2004[ ). Both chiral and 
continuum extrapolations have been taken, though the pion masses in the simulation 
are rather heavy. 

The / = 2 TTTT scattering length has been calculated in the 2-|-l-flavor mixed action 
lattice QCD, using the domain- wall valence quarks with the asqtad-improved staggered 
sea quarks for mf° ~ 294, 348 and 484 MeV at a ~ 0.125 fm ^eane et al., 2006[ 



Beane et al, 2008). The scattering phase shift has also been calculated at ~ 544 
MeV and ~ 484 MeV. 

Recently the / = 2 tttt scattering length has been calculated in 2-flavor twisted 
mass lattice QCD for pion masses ranging from 270 MeV to 485 MeV at a ~ 0.086 
fm( |Feng et al., 2010[ ). The lattice spacing error is estimated at a ~ 0.067 fm for one 
pion mass. 

The P-wave scattering phase shift for the / = 1 tttt system has been calculated in 
2-flavor lattice QCD with an improved Wilson fermion at a = 0.22 fm in the laboratory 
system( pVoki et al., 2007[ ). Since mT^/mp ~ 0.41 in this calculation, the decay width of 
p meson can be estimated from the scattering phase shift. 



3 

Nuclear Potential from Lattice QCD 



In the previous section, we have explained the Liischer's method to extract the scat- 
tering phase shift from the two particle energy in the finite box, considering the tttt 
case as an example. To show the relation between the phase shift and the two particle 
energy in the finite box, we use the fact that the BS wave function in the large sepa- 
ration such that r > R, where R is the interaction range between two particles in the 
infinite volume, satisfies the free Schrodinger equation (the Helmholtz equation) with 
the periodic boundary condition. 

In this section, instead of the large distance behaviour, we consider the short dis- 
tance properties of the BS wave function, from which we define the " potential" between 
two particles. We mainly consider the NN potential, though the method in this section 
can be applied to any two particles in principle. 

3.1 Strategy to extract potentials in quantum field theories 

In this subsection, we describe the strategy to extract the NN potentials in QCD pshii et ai, 2007[ 



Aoki et ai, 2008 Aoki et al., 20106) 



As a preparation, we introduce the T-matrix of the NN scattering below the NNtt 
inelastic threshold. The 4x4 T-matrix component for a given total angular momentum 
J is decomposed into two 1x1 submatrices and one 2x2 submatrix as( |Ishizuka, 2009[ 
|Aoki et ai, 20106D 

1,7 _ I n t7 



T'^\ T^tj^^^, ^0ix2 I (3.1) 

02x1 02x1 T(Lj^i,s=l. 



where / is the orbital angular momentum between two nucleons and s is the total spin. 
The unitarity tells us that 



Tlj, = fj,, TiLj^,„,^, = 0{k) ( ^'^^i^i ^ ) 0~\k) (3.2) 



with 



k \^sine,7(fc) cose,/(fc) J' ^ ' 

where 5is{k) is the scattering phase shift, whereas ej(fc) is the mixing angle between 
/ = J ± 1. Here the total energy of the two nucleons is given by 2Ek = 2^Jk'^ + in 
the center of mass frame. 
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Let us start describing the strategy to extract the potential in QCD. We first define 
the BS amplitude for two nucleons in the center of mass frame as 

(^f^(r) = (O|T{iV„(y,O)iV0(x,O)}|k,Sa,-k,s,;in), (3.4) 

where the relative coordinate denoted as r = x — y, the special momentum and the 
helicity for incoming nucleon is denoted by (k, Sa) or (— k, Sb), the total energy E ~ 2Ek 
with k = |k|. The local composite nucleon operator is given by 

Nlix) = 6'^''^gS'-^(x)g^^^(z)(zr2),.(C75)^^gf (x) 

= e^'^q^^fix) [q^(cc)zT2C75g^(x)] , (3.5) 

where q^'-^ is a quark field with the color index a, the fiavor index / and the spinor index 
a. Here repeated index assumes a sum, C = 7274 is the charge conjugation matrix 
and iT2 acts on the flavor index. Unless necessary, the flavor indices are implicit. 
Let us briefly consider the meaning of the BS wave function. By writing 



T{iV(y,0)iV(x,0)} = V I ^{2N,n{NN),mn,E\U^ir,E) 



^ J 2E 

n.rn— 



dE 



(3.6) 



where \2N , n{N N) , mn , E) is an in-state containing two nucleons, n pairs of nuclcon- 
antinucleon and m pious with the total energy E, we see that ^p^{v) = /oo(r, E). (Our 
normalization is (27V, n(7V7V), m7r, i;| |2iV, n'(7V7V), mV, £")) = 2E5{E - E')5nn-5nim' ■ 
) Therefore the BS wave function ^p^{v) is an amplitude to find the in-state \2N,E) 
inr{7V(y,0)7V(x,0)}|0). 

As in the case of tttt, the asymptotic behaviours of the BS wave function at r = 
|r| > i?, where R is the interaction range of two nucleons, agree with those of the 
scattering wave of the quantum mechanics! Ishizuka, 2009 Aoki et ai, 2010&). The BS 
wave function for a given total angular momentum J has 4 components. For example, 
the BS wave function for Z = J, s = at large r becomes 

(r; k) ZYjj^{nr)e'^-"''-''^ (jj(fcr) cosdjo{k) + nj{kr) sindjo{k)) 

r>R 

piSjoik) 

~ ZYjj,{n,)— sin(fcr + Sjo{k) - nJ/2). (3.7) 

kr 

where is the z component of the total angular momentum. See Refs. Pshizuka, 2009[ 
|Aoki et al, 2010"6| for more details. This shows that the BS wave function can be 
regarded as the NN scattering wave. 



Now we define the non-local NN "potential" through (r)(Ishii et ai, 2007 



Aoki et ai, 2008 Aoki et al., 201061 as 



{'k " ^"''^''^ " / "^'^ Uap,^s{^,y)v^s{y): Ho = (3.8) 

where ^ = 111^/2 is the reduced mass of the two nucleons. It is noted that C/(x, y) 
is non-local but energy-independent and this potential is equivalent to the local but 
energy dependent potential F(r, k), which is defined by 
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K{r, k) = (^^ - Ho^ ^{r, k) = F(r, k)ipi^, k), (3.9) 

where we write (/^(r, k) = ip^{r). To see the equivalence, we construct the dual basis 
(p(k, r) as 

(^(k,r) = J d^p 77-i(k,p)<^(r,p), (3.10) 
where the metric is given by 

7?(k,p) = J dV<^(r,k)(^(r,p). (3.11) 
It is easy to see that the dual basis satisfies 

J d^y ^(k,r)^(r,p) = J(3)(j^_p) (3^^2) 
d3p^(x,p)^(p,y) = (3^13) 

Using the dual basis, we obtain the non-local potential form the local one as 

C/(x, y) = y d^p A'(x, p)(^(p, y)= J d^p ^(x, p)^(x, p)(^(p, y). (3.14) 

This establishes one to one correspondence between the non-local but energy-independent 
potential [/(x,y) and the local but energy-dependent one y(r, k). 

The equivalence also tells us that we need to know the BS wave function at all 
energies to completely construct J7(x, y). Although this is principle possible, it is in 
practice very difficult. We therefore consider the following derivative expansion of 
C/(x,y). 

(7(x,y) =y(x,V)5(3)(^_y)^ (3 15) 

The structure of the V{r, V) can be determined as follows( Okubo and Marshak, 1958] ). 
The most general (non-rclativistic) potential is parameterized as 

y(ri,r2,Pl,P2,CTi,CT2,Ti,f2,0 (3.16) 

where r^, p^. (Ti and fi arc the coordinate, the momentum, the spin and the isospin of 
the i-th nucleon, respectively and the t is the time. There are several conditions this 
potential should satisfy. 

1. Probability conservation implies the hermiticity of the potential, = V . 

2. Encregy conservation imply the t independence while the momentum conservation 
says that the potential depends on the combination r = ri — r2 only. 

3. Galilei invariance tells us that the potential contains p = Pi — P2 only. From these 
three conditions, we have V = y(r, p, 172, fi, f^). 



Strategy to extract potentials in quantum field theories 21 



4. The total angular momentum conservation implies that V is invariant under J = 
L + S with the orbital angular momentum L = r x p and the total spin S = 

(^1 + CT2)/2. 

5. The potential should be invariant under parity, (r, p, cfi) — > (— r, — p, ct,). 

6. The potential is invariant under time-reversal, (r, p, ai) — (r, — p, — CTj). 

7. Quantum statistics of the exchange of two nucleons implies the invariance of the 
potential under (r, p, cti, 1T2, ti, T2) ^ (^r, — p, (72, (Ti, T2, 

8. From isospin invariance, V contains only 1 • 1 or ri • T2 in the isospin space. 

9. The potential has only a'^a^ terms with (n, m) = (0, 0), (1, 0), (0, 1), (1, 1). The 
other higher order terms can be reduced to these terms because of the property 
that crV^' = +ie'i^(j^. 

The terms which contain Pauli matrices and satisfy the above conditions are con- 
structed as 

CTi • CT2, {Si + CT2) • {Si ■ r)(CT2 • r), {Si ■ p)(ct2 • p), {Si ■ L)(ct2 • L), (3.17) 
which are customarily reorganized as 

<7l-(?2, S'12 = 3((7i • f)((72 • f) - (7l • (T2, L-S, 

Pi2 = (ai •p)(a2 -P), Wi2 = Qi2-\Si-S2p, (3.18) 

where 512 is called the tensor operator, and 

1 



/12= 2 



^ {Si ■ L){S2 ■ L) + {S2 ■ L){Si • L),J , (3.19) 
Finally we obtain 

V=Y, V\v,p,Si,S2)PJ (3.20) 



/=1,2 

where 

V' = + Vj{Si ■ S2) + Vls{L ■ S) + \{Vi, S12} + \{V^, P12} + \{V^,Wi2} 

(3.21) 

with coefficient functions = (r^, p^, L^) for / = 0, 1 and X = 0, a, T, LS, P, W. 
Here PJ is the projection operator to the state with the total isospin /, given by 

P[=, ^\-ri-f2, PU = ^ + Ti • f2 . (3.22) 

The anticommutators in Eq. p.2ip are necessary to make the potential hermitian, since 
'5'i2, P12 and W12 do not commute with the scalar potentials V^(r^, p^, L^). 
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From the general consideration above, we have, for example at 0(V), 



F(r, V) = ^ \vi{r) + y/(r)ai • + Vi{r)S,2 + Vls{r)L ■ s] PJ 



+ 0(V2). 



(3.23) 



This form of the potential has often been used in nuclear physics. Note that the first 
3 terms are 0(1) while the LS potential is 0(V). 

This is the strategy to define and extract the NN potential in QCD. There are 
two important and mutually-related remarks. 

The leading order potential in the derivative expansion is nothing but the local 
potential l/(r,k). By construction, this (local) potential reproduces the correct phase 
shift 5{k) at k = |k|, while it is not guaranteed that this potential gives the correct 
phase shift at different k' = |k'|(7^ k). This means that the local potential V{r,'k) and 
V^(r,k') may differ, and this energy dependence of the local potential gives a measure 
for the non-locality of the non-local but energy independent potential U (x, y) because 
of the equivalence between V and U. If the the first order in the derivative expansion 
is good at low energy, we expect that energy dependence of the potential V is expected 
to be small at small k. 

Secondly it should be mentioned that the potential U defined through the BS 
amplitude of course depends on the choice of the interpolating field operators N{x). 
In principle, one may choose any (local) composite operators with the same quantum 
numbers as the nuclcon to define the BS wave function. Different choices for the 
nuclcon operator give different BS wave functions, which may leads to different NN 
potentials, though they all gives the same scattering phase shift. While the potential 
is not an physical observable in this sense, it does not mean that it is useless, however. 
The strategy in this report gives one specific scheme for the NN potential in QCD, 
which is defined through the BS amplitude constructed from the local nucleon field 
without derivatives. This is quite analogous to the situation for the running coupling in 
QCD. Although the running coupling is scheme-dependent, it is useful to understand 
and describe the deep inelastic proton-electron scattering. Let us make this analogy 
more concrete. A physical observables is the scattering data in both cases, the deep 
inelastic scattering or the NN scattering. An example of the physical interpretation 
is the almost free partons in proton for the case of the deep inelastic scattering, while 
it is an existence of the repulsive core for the case of NN scattering. An theoretical 
explanation for the phenomena is the asymptotic freedom of the QCD running coupling 
for the free partons, while no valid theoretical explanation exists so far for the repulsive 
core. In this report, we introduce one definite scheme for the potential based on QCD, 
in order to show an existence of the repulsive core. Although the choice of the scheme 
is irrelevant in principle, it is better to use a "good" scheme in practice. In the case 
of the running coupling, good convergence of the pertubative expansion may give one 
criterion, though the popularly used MS coupling may not be the best one for this 
criterion. In the case of the NN potential, on the other hand, good convergence of 
the derivative expansion may give a criteria for a "good" potential. In other words, 
the good potential is almost local and energy-independent. The NN potential which 
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is completely local and energy-independent at all energy range is therefore the best 
one. It is also unique if the inverse scattering method holds for the A'' A'' case. 

3.2 Extraction of the BS wave function on the lattice 

In this subsection, wc explain how to extract the BS wave function from correlation 
functions on the lattice. For simplicity, we here consider the I = state, namely the 
S'-state. 

The BS wave function on the lattice with the lattice spacing a and the spatial 
lattice volume is extracted from the 4-point correlation function, by inserting the 
complete set of the QCD cigenstates in the finite box as, 

Ga;3(x,y,t-io; J^) = (0|n^(y,tK(x,i)Jp„(to; J^)|0) (3.24) 

OO 

= ^ A„(0|n^(y,tK(x,t)|£;„)e-^"(*-*°) (3.25) 

ji=0 

> Aaipap{r;J^), r = x-y (3.26) 

with the matrix element An = {En\Jpn{0)\0) , where p (n) is the proton (neutron) 
interpolating operator, p = iV" ( n = -/V*), \En) is the QCD eigenstate with the 
baryon number 2 and the total energy En = 2^k\ + m^. The state created by the 
source Jp„ have the conserved quantum numbers, {J,Jz) (total angular momentum 
and its z-component), /(total isospin) and P (parity). To study the A'' A'' potential in 
the JP = 0+ with / = 1 (iSo) channel and the = 1+ with 7 = (^Si and ^Di) 
channel, a wall-source located at t = to with the Coulomb gauge fixing only at t = to 
is used, 

Jpn{t,, J^) = P^^ [pr"(io)n^^"(to)] (3.27) 

where p^'^^^ito) and n^^''(to) are obtained by replacing the local quark fields q{x) in 
N{x) by the wall quark fields, 

9™"(io)=5^9(x,<o). (3.28) 

X 

By construction, the source operator ea. (|3.27p has zero orbital angular momentum 
at t ~ to, so that states with fixed {J,Jz) are obtained by the spin projection with 
(s, Sz) = (J, Jz), e.g. Pp~° = (cr^)/3a and -P|~^'''''~° = {o'^)/3a- Note that the I and s are 
not separately conserved, so that the state created by the source Jpn{to] l"*") becomes 
a mixture of the / = (S-state) and I = 2 (D-state) at later time t > to- 

The BS wave function in the orbital S state is then defined with the projection 
operator for the cubic group P^ with the irreducible representation R and that for 
the spin P'^ as 

^(r;iSo)=P^^F'^=Vr;0+)^ 2i E Pap\M9-'r;0+) (3.29) 

S60(3,Z) 
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'Si)=P^^P«=V(r;l+)^^ ^ P„7VaM3-'r;l+), (3.30) 



geO(3,Z) 



where the summation over g € 0(3, Z) is taken for the cubic transformation group 
with 24 elements to project out the I = component in the Af representation, and 
contributions from the higher orbital waves with I > 4 contained in the rep. are 
expected to be negligible at low energy. 

Form these BS wave functions, we can construct the local potentials at the leading 
order of the derivative expansion. For the ^Sq channel, the central potential becomes 

Voir;' So) ^ Vt'{r) + Vj=\r) ^ ^ + — ^'f^^''^ (3.31) 

niN rriN ip{r;'- bo) 

while for the "^Si channel there are two independent terms, Vc{r,^ Si) = Vj/^"(r) — 
3V/^°(r) and V^^^{r), at the leading order. For a while we ignore Vp and define the 
effective central potential as 

Vf(r:^S.)^^ + ^^^^. (3,32) 

mN rriN ip{rfbi) 

where the "effective" central potential means that it includes the effect of the tensor 
potential Vr as the second order pertubation. 

It is noted here that k'^ is determined from the total energy Eq of two nucleons as 

3.3 Tensor potential 

While the central potential acts separately on the S and D components, the tensor 
potential provides a coupling between these two components. We therefore consider a 
coupled-channel Schodinger equation in the = 1+ channel, in which the BS wave 
function has both S-wave and D-wave components as 

{Ho + Vcir; 1+) + VT{r)Si2) <^(r; 1+) = —^{r; 1+). (3.33) 

rriN 

The projections to the S-wave and D-wave components similar to eq. p.30p are defined 

by 

P^„^(r) = P^^^„^(r;l+), (3.34) 
QiPc^pir) = (1 - P<)^«^(r; 1+). (3.35) 

Here both Pipa/i and Qipap contain additional components with I > 4 but they are 
expected to be small at low nenergy. 

Multiplying P and Q to eq. (|3.33p from the left and using the properties that Hq, 
Vc{r; l"*") and Vrir) commute with P and Q, we obtain 

Ho[P^]ir) + Vcir; l+)[P¥>](r) + VTir)[PSi2'p]ir) = — [P(^](r) (3.36) 

TTlN 
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Fig. 3.1 The A''A'^ wave function in ^So and '^Si cliannels at = 529 MeV, measured at 
(t — to) = 6a. Tlie inset is a tliree-dimensional plot of tiie wave function ip{x,y,z = 0;^ So). 

Ho[Q^]{r) + Voir; l+)[Q^]{r) + VT{r)[QS,2ip]{r) = —[Q^Kr) (3.37) 

niN 

where, for simplicity, the spinor indices, a and (3 are suppressed. Note that 15*12 does 
not commute with P or Q. 

By solving these equations for (a, (3) = (2, 1) component, we finally extract 

Voir; 1+) = — - -1- {[QSi2^]2i{r)Ho[P^]2iir) - [PSi2'p]2i{r)Ho[Q^]2iir)} 



Vrir) 



rriN D{r) 
1 



D{r) 



{[Q^]2i{r)H^[P^]2i{T) - [P^]2i(r)i/o[Q^]2i(r)} 



where 

D{r) = [P<p\2i{r)H^[QSi2V]2i{r) ~ [Q(^]2i (r)iJo[P5i2(^]2i(r). 
Note that the efi'ective central potential is expressed as 



rriN [Pip\2i{r) 



(3.38) 
(3.39) 

(3.40) 

(3.41) 



with iJo = —'V'^/mpf. 

3.4 Results in lattice QCD 

The first result for the NN potential in lattice QCD based on the strategy in the 
previous subsections appeared in Ref.( |Ishii et ai, 20Q7[ ), where the (effective) central 
potential has been calculated for ^Sq ('^Si) channel in quenched QCD simulations at 
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Fig. 3.2 The central potential in the ^So channel and the effective central potential in the 
^Si channels at = 529 MeV. 

the lattice spacing a ~ 0.137 fm and the spacial extension L ~ 4.4 fm. More details of 
numerical simulations can be found in Ref. ( jlshii et al., 2007] ). 

Fig l3. II shows the BS wave function in -'^So and '^Si channels at TOtt = 529 MeV and 
fc^ ~ 0. which is measured at i — = 6a. The wave functions are normalized to be 1 
at the largest spatial point r = 2.192 fm. 

The reconstructed central and effective central potentials in the "'^Sq and "^Si chan- 
nels at TOtt = 529 MeV from the BS wave functions with the formulae p.3ip and (|3.32p 
are shown in Fig |3.2l The overall structure of the potentials are similar to the known 
phenomenological NN potentials discussed in SecHJ namely the repulsive core at short 
distance surrounded by the attractive well at medium and long distances. The figure 
also shows that the interaction between two nucleons is well switched off for r > 1.5 
fm, so that the condition R < L/2 ~ 2.2 is satisfied. 

To check the stability of these potentials against the time-slice adopted to define 
the BS wave function, the t-dependence of the ^Sq potential for several different values 
of r is shown in Fig |3.3l at tUt^ = 529 MeV. In this case, the choice t — to = 6a for 
the extraction of Vc(r) is large enough to assure the stability within statistical errors, 
which indicates the ground state dominance at this t. 

The NN potentials in the ^Sq channel are compared among three different quark 
masses in Fig l3.4l As the quark mass decreases, the repulsive core at short distance 
and the attractive well at medium distance becomes stronger simultaneously. 

The '^Si and ^Di components of the BS wave functions obtained from = 1+, 
Jz = state at 171^^ ~ 529 MeV and A:^ ~ are plotted in Fig l3.5r a). according to 
eqs. p.34p and p.35|) . Note that the '^Di wave function becomes multi-valued as a 
function of r due to the its angular dependence. It is expected that (a, /?) = (2, 1) spin 
component of the D-statc wave function for = 1+ and = is proportional to the 
^20(0, </>) oc Scos'^ 6' — 1. As shown in Fig |3.5f b). the D-state wave function, divided by 
the 120(^7 "Z"); becomes almost single-valued, so that the D-wave component is indeed 
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Fig. 3.3 The t dependence of the potential at r = 0, 0.14, 0.19, 1.37, 2.19, 0.69 fm from top 
to bottom in the ^So channel at uitt — 529 MeV. 
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Fig. 3.4 The central potentials in the ^So channel at three different quark masses. 

dominant in Qip{r). The central potential Vc{r;l~^) and the tensor potential Vrir) 
together with the effective central potential V^^{r;^ Si) in the 3Si channel are plotted 
in Fig l3.6l Note that VQ^{r;^ Si) contains the effect of Vrir) implicitly as higher order 
effects through the process such as ^Si ^3 ^3 ^^le real world, Vf{r;^ Si) 

is expected to acquire large attraction from the tensor force, which is reason why 
the bound state for a deuteron exist in the '^5'i while no bound states appears for a 
dineutron. As seen from Fig |3.61 the difference between Vc{r; 1+) and V^^^rf Si) is 
still small in this quenched simulations due to relatively large quark masses. 

The tensor potential Vr (r) in Fig l3.6l is negative for the whole range of r within 
statistical errors and has a minimum at short distance around 0.4 fm. If the tensor 
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Fig. 3.6 The central potential Vc{r;V^) and the tensor potential Vrir) obtained from the 
J+ BS wave function at = 529 MeV, together with Vc^{r-^ Si). 

force receives significant contribution from the one-pion exchange as expected from 
the meson theory, Vt would be rather sensitive to the change of the quark mass. As 
shown in Fig l3.7i indeed the attraction of Vrir) substantially increases as the quark 
mass decreases. 

At present potentials are determined at leading order of the derivative expansion, 
and examples presented so far are extracted from lattice data taken at fc ~ 0. If the 
higher order terms such as VLs{f)L ■ S becomes important, the LO local potentials 
determined at fc > are expected to be different from the one at fc ~ 0. From such k 
dependence of the LO local potentials, some of the higher order terms can in principle 
be determined. A lattice QCD analysis on the k dependence has been recently carried 
out by changing the spatial boundary condition of the quark field from the periodic 
to the anti-periodic ones, which corresponds to the change from fc ~ MeV to fc = 
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Fig. 3.7 Quark mass dependence of the tensor potential Vrir). 

y'3(7r/L)^ ~ 250 MeV. In Fig l3.8l the local potential for the ^5*0 channel obtained at 
k ~ 250 MeV is compared with the one at fc ~ in quenched QCD at a = 0.137 fm 
and ttItt = 529 MeV. As seen from the figure, the k dependence of the local potential 
turns out to be very small for every r within statistical errors. Namely the non-locality 
of the potential with the choice of the local interpolating operator for the nucleon is 
small, and the present local potential at the LO can be used to well describe phsyical 
obscrvables such as the phase shift (So(fc) from fc ~ to fc ~ 250 MeV without significant 
modification, at least in quenched QCD at a = 0.137 fm and TOtt = 529 MeV. This 
also indicates that the definition for the potential through the BS wave function with 
the local nucleon operator is a " good scheme" . 



30 Nuclear Potential from Lattice QCD 




Fig. 3.8 A comparison of central potentials at fc ~ (PBC, blue) and at fc ~ 250 MeV 
(APBC, red) for the ^So state in quenched QCD at a = 0.137 fm and = 529 MeV. 
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Repulsive core and operator product 
expansion in QCD 



As shown in the previous section, the lattice QCD calculations shows that the A^A^ 
potential defined through the BS wave function has not only the attraction at medium 
to long distance, which has long well been understood in terms of pion and other heav- 
ier meson exchanges, but also a characteristic repulsive core at short distance, whose 
origin is still theoretically unclear. A recent attempt( |Aoki et ai, 2010"a| to theoreti- 
cally understand the short distance behavior of the potential in terms of the operator 
product expansion(OPE) is explained in this section. 

4.1 Basic idea 

Let us first explain the basic idea. We consider the equal time BS wave function defined 

by 

V^f ^(r) = (0|O^(r/2, 0)OB(-r/2, 0)\E) (4.1) 

where \E) is some eigen-state of a certain system with the total energy E, and Oa, 
Ob are some operators of this system. We here suppress other quantum number of the 
state \E) for simplicity. The OPE reads 

0A(r/2, 0)OB(-r/2, 0) ^ ^ D%ir)Oc{0, 0), (4.2) 

c 

from which we have 

^ABir) c^J2DAB{rmOc{0,0m. (4.3) 

c 

It is noted that r dependence appears solely in D'^g (r) while the E dependence exists 
only in (0|Oc(0, 0)|£^). Suppose that the coefficient function of the OPE behaves in 
the small r = |r| limit as 

D'^Bir) ^ r^^{~\ogrf-fc{e,ct>) (4.4) 
where 0, are angles in the polar coordinate of r, the BS wave function becomes 

^fB(r) ~ ^r"^(-logr)^^/c(0,0)i^c(i?), Dc{E) = {0\Oc{0,0)\E) (4.5) 
c 
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in this limit. The potential at short distances can be known from this expression. For 
example, in the case of the Ising field theory in 2-dimensions, the OPE for the spin 
field (T is given by 

cr(a;,0)cr(0,0) ~G(r)l+cr3/4C)i(0) + ..., r = (4.6) 

where Oi (x) (=: 'ip^{x) : in terms of the free fermion fields) is an operator of dimension 
1, which leads to 

ip{r,E)'^r^'^D{E)+0{r'^''^), D{E) = c{{)\Oi{fd)\E) , (4.7) 



where \E) is a two-particle state with energy E = 2\/fc2 + m^. From this expression 
the potential becomes 

= + (4.8) 

mtp{r, E) 16 mr^ 

in the r limit. The OPE predicts not only the r^"^ behavior of the potential at short 
distance but also its coefficient —3/16. Furthermore the potential at short distance does 
not depend on the energy of the state in this example( |Aoki et al., 20096] ). 

As will be seen later, the dominant terms at short distance have ac = 0. Among 
these terms, we assume that C has the largest contribution such that Pc > Pc for 
VC" ^ C. Since, as will be also seen later, such dominant operators with ac = 
mainly couples to the zero angular momentum (L — 0) state, let us consider the BS 
wave function with L ~ 0, given by 

Vfsir) {-logrf^ Dc{E) + V {~\ogrfc'Dc'{E). (4.9) 



Using 



(-logr)''-! r 13 ~ I 



1 - 



(4.10) 



no ^-1^7^^, (4.11) 



lo 

we obtain the following classification of the short distance behavior of the potential. 

1. /3c 7^ 0: The potential at short distance is energy independent and becomes 

r'^{— logr) 

which is attractive for /J^ > and repulsive for /3c < 0. 

2. ;9c = 0: In this case the potential becomes 

where /3c' < is the second largest one. The sign of the potential at short distance 
depends on the sign of Dc'{E) / Dc{E). 

On the lattice, we do not expect divergence at r = due to lattice artifacts at 
short distance. The above classification hold at a <C r <C 1/Aqcd, while the potential 
becomes finite even at r = on the lattice. 



Renormalization group analysis and operator product expansion 33 



4.2 Renormalization group analysis and operator product 
expansion 



In A^/-flavor QCC regularized in _D = 4 — 2e, bare local composite operators O^' (x) 
are renormalized as 

0^r\^) = ZAB{g,e)0^°\x), (4.13) 

We here ignore the contribution from quark mass terms, which generates less singular 
terms in the OPE. Throughout this section, summation of repeated indices is assumed. 
The meaning of this formula is that finite results arc obtained if we insert the right 
hand side into any correlation function, provided the QCD coupling and the quark and 
giuon fields are appropriately renormalized. We first consider an n-quark correlation 
function without operator insertion Gn\goj^) (Here the dependence on the quark 
momenta and other quantum numbers suppressed.), which is renormalized as 

G(-")(5,A^) = Z-"/^(g,e)Gio)(5o,6), (4.14) 

where the coupling renormalization is given by 

g'„^^L''Z,{g^e)g\ (4.15) 

The renromalization constant Zi in the minimal subtraction (MS) scheme has only 
pure poles terms as 

^l(<7^e) + ^ + o(/) (4.16) 

where 

Similarly the quark field renormalization constant is given by 

Z^^(5,e) = l-^ + 0(/). (4.18) 



The gluon field renormalization constant is also similar but is not necessary for our 

3) / 
.■.A\ 



purpose. Similarly an n-quark correlation function with operator insertion G^^}j^{go, e) 



is renormalized as 

G<-:^^\g, /i) = ZMg, Z^^'ig, ^) G%{g,. e) , (4.19) 

where 

(1) 2 

ZABig, e) = Sab - + 0{g^). (4.20) 

The renormalization group (RG) equations are obtained from the simple fact that 
bare quantities are independent of the renormalization scale /i. Introducing the RG 
differential operator 
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(4.21) 

the RG equation for n-quark correlation functions can be written as 

+ ^7F(g)} G(:™'(5,m) = 0, (4.22) 

where the RG beta function is 

pD{g) ^ = -—i^T^ = -eg - Pog' - Pig' + 0{g'), (4.23) 

^ + 2 

while the RG gamma function for quark fields is 

Mg) - Mg, e) = 7F0 <?' + 0(/). (4.24) 

Note that Poig) differs from /34(g) only by —eg, and therefore has the smooth limit 
to = 4. The RG invariant A parameter satisfies D{^)A = with the Ansatz that 

A = fie^^o\ (4.25) 

The solution is the lambda-parameter in the MS scheme, Ams, if the arbitrary inte- 
gration constant is fixed by requiring that for small coupling 

/(.) = -^-|.ln(/3o.^)+0(,^). (4.26) 

Finally the RG equations for n-quark correlation functions with operator insertion are 
of the form 

{Di^,) + ^7F(.g)} Gi;A^(5,A^) -7AB(g)GS^(.9,/.) = 0, (4.27) 

where 

lABig) = -ZAcPnig, e)^ = + 0(/). (4.28) 

Let us consider the OPE 

Oi(y/2)02(-2//2) DB{y) Os(0). (4.29) 

where Oi and O2 are nucleon operators and the set of operators Ob are local 6-quark 
operators of canonical dimension 9 and higher. All operators in (|4.29p are renormalized 
ones, but from now on we suppress the labels As we will see, the nucleon operators 
are renormalized diagonally as 

0i^Zi{g,e)0f\ 02-^2(.9,e)0f , (4.30) 

and the corresponding RG gamma functions are defined by 
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7i.2(.9) - Mg, e) = ^['jg' + 0(/). (4.31) 

Comparing (|4.29p with its bare version, 

0['\y/2)0i"\~y/2) ^ 7^^"' (y) O^"' (0), (4.32) 
we can read off the renormahzation of the coefficient functions as 

DB{y) ^ Zi{g,e)Z2{g,e)D^°\y) Z^},{g,e) (4.33) 
so the the RG equation becomes 

D{ii)DBig, ti,y) + DAig, ^i,y)lAB{g) ^ 0, (4.34) 
where the effective gamma function matrix is defined as 

lAsig) = lAB{g) - [71(5) +72 (,9)] 5ab- (4.35) 

Hereafter we assume the dimensionless coefficient functions, which can be written as 
DAig, fJ-, y) — DA{g] A^) with r = \y\. For the case of operators with higher canonical 
dimension 9 + a the coefficients are of the form r" times dimensionless functions 
and the analysis is completely analogous and can be done independently, since in the 
massless theory operators of different dimension do not mix. In the full theory quark 
mass terms are also present, but they correspond to higher powers in r and therefore 
can be neglected. 

To solve the vector partial equation (|4.34p . we introduce UAB{g), the solution of 
the matrix ordinary differential equation 

Pig) ^ UABig) = lAcig) UcB{g) (4.36) 

and its matrix inverse UABig)- With this solution, DBig'^^J-r) can be easily obtained 
as 

Deig; fir) = ^^(Ar) UABig), (4-37) 

where the vector Fa is RG-invariant. Now the running coupling g is introduced as the 
solution of the equation 

fig) = fig) + Hf^r) = In(Ar) . (4.38) 

Note that ^ is a function of r but does not depend on fi. Since Fb is RG invariant, we 
can evaluate it at /i — 1/r as 

FBiAr) = DAig;l)UABig)- (4.39) 

By definition g = g at — 1/r. Since, because of asymptotic freedom (AF), for ?' — > 
also g — as 

Fb can be calculated perturbatively. 
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Putting everything together, the operator product expansion (|4.29p can be rewrit- 
ten as 

Oi(2//2)02(-y/2) ~ FeiAr) Ob(0), (4.41) 

where 

dB^UBc{9)0c. (4.42) 

There is a factorization of the operator product into perturbative and non-perturbative 
quantities: FsiAr) is perturbative and calculable (for r — > 0) thanks to AF, whereas 
the matrix elements of Ob are non-perturbative but r-independent. 
The coefficient functions have the perturbative expression 

DA{g; fir) = Da-m + g^DAaif^r) + 0(/), (4.43) 

and the basis of operators can be chosen such that the 1-loop mixing matrix is diagonal: 

jABig) = 2/3o Pa 5' Sab + 0(/). (4.44) 

In such a basis the solution of ()4.36p in perturbation theory takes the form 

UABig) = {Sab + RABig)} g-^^", (4.45) 

where RAB{g) — 0{g^), with possible multiplicative \ogg^ factors, depending on the 
details of the spectrum of 1-loop eigenvalues Pa- An operator Ob first occurring at 
^B-loop order on the right hand side of ()4.29p has coefficient FB(Ar) with leading short 
distance behavior 

FBiAr) « i?s/,(l)5'^'""^"^ « DbjJI) (-2/3o In(Ar))^--'- . (4.46) 

In principle, an operator with very large /3b, even if it is not present in the expansion 
at tree level yet, might be important at short distances. This is why it is necessary to 
calculate the full 1-loop spectrum of Pb eigenvalues. As we will see, no such operators 
exist for the two nucleon case, and therefore operators with no n- vanishing tree level 
coefficients are dominating at short distances. The corresponding coefficient functions 
have leading short distance behavior given by 

FB{Ar) « Db,o (-2/3o HMf" ■ (4-47) 

4.3 OPE and Anomalous dimensions for two nucleons 

The general form of a gauge invariant local 3-quark operator is given by 

(x) ^ Bif^ix) ^ e^^'^ql:f{x)q''^%x)q^^\x) , (4.48) 

The color index runs from 1 to = 3, the spinor index from 1 to 4, and the flavor 
index from 1 to Nf. A summation over a repeated index is assumed, unless otherwise 
stated. Note that B^J^^^ is symmetric under any interchange of pairs of indices (e.g. 

^Q^7 ~ ^^Q-y) because the quark fields anticommute. For simplicity wc sometimes use 
the notation such as = fgh and F = a/37 as indicated in (|4.48p . 
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The nuclcon operator is constructed from the above operators as 

Blix) = (P+4)„„, B^?^^(C75)0.(^r2)^" , (4.49) 

where P+4 = (1 + 74)72 is the projection to the large spinor component, C = 7274 
is the charge conjugation matrix, and T2 is the Pauh matrix in the flavor space (for 
Nf = 2) given by {1x2)-^^ — s^^ . Both C75 and 1x2 are anti-symmetric under the 
interchange of two indices, so that the nucleon operator has spin 1/2 and isospin 1/2. 
Ahhough the exphcit form of the 7 matrices is unnecessary in principle, we find it 
convenient to use a chiral convention given by 

o''^ ' 74 = (^^ , 75 = 71727374 = J ■ (4-50) 

As discussed in the previous subsection, the OPE at the tree level (generically) 
dominates at short distance. The OPE of two nucleon operators given above at tree 
level becomes 

B^x + y/2)Bl{x yl2) = i3£(x)i?^(x) + ^ {a^[i?/(x)]i?|(x) - Bi{x)d,[Bl{x)]] 

+ O(y^). (4.51) 

For the two-nuclcon operator with cither the combination [a/3], {fg} {S — 0) or the 
combination {a/?}, [fg] {S = 1), terms odd in y vanish in the above OPE, so that only 
even L contributions appear. These 6-quark operators are anti-symmetric under the 
exchange (a, /) O (/3,5). On the other hand, for two other operators with ([a/3], [fg]) 
or ({a/3}, {/.?}), which are symmetric under the exchange, terms even in y vanish in 
the OPE and only odd L's contribute. 

Knowing the anomalous dimensions of the 6-quark operators appearing in the 
OPE, which will be calculated later in this subsection, the OPE at short distance 
(r = |y| <c 1, y4 = 0) becomes 

Biix + y/2)Blix~y/2) ^ ^ c^(r)0^«/(x) + ^ ds(r)/0^^^;f (x) + 0(r2) , 

A B 

(4.52) 

where the coefRcient functions behave as 

CA{r) ~ (- logr)'3-^ , deir) ~ (- logr)'^^ , (4.53) 

and (3a, B are related to the anomalous dimensions of the 6-quark operators O^f^'^ and 

of those with one derivative O;^^'^. The wave function defined through the eigenstate 
\E) is given by 

^r'^(y) = {0\Biix + y/2)B$ix - y/2)\E) ^ ^ CAirMOi^^ {x)\E) (4.54) 

A 

for the anti-symmetric states, while 
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V°j,^'iy) = {0\Biix + y/2)Bf,{x-y/2)\E) ^Y.Mr)y'{0\Oi$^!',ix)\E) (4.55) 

B 

for the symmetric states. Hereafter we consider only 6-quark operators without deriva- 
tives and calculate the corresponding anomalous dimensions. 

The renormalization factor Zx of a fc-quark operator X = [g*^] is defined through 
the relation 

[q'^Y^^^Zx[q'o]^ZxZ'J\''], (4.56) 

where qoiq) is the bare (renormalized) quark field. The wave function renormalization 
factor for the quark field is given at 1-loop by 

Z^ = l+.g^zW, 4^^ = -^ (4.57) 

where A is the gauge parameter and Cp ~ ^2N^ ■ 

At 1-loop the renormalization of simple fc-quark operators without gauge fields is 
given by the divergent parts of diagrams involving exchange of a gluon between any 
pair of quark fields. The 1-loop correction to the insertion of an operator q^f {x)q'p^ (x) 
in any correlation function involving external quarks is expressed as the contraction 
of 

CHx)ql'''{x)^_ [ d^V d^z A^^{y)A^{z)[qf^ {y)^gT^l,q^' {y)W {z)^gT^l.q'' [z)] 



(4.58) 



where tr T^^T^ 5"^^ 12 in our normalization. Since two identical contributions cancel 
the 2! in the denominator, the contraction at 1-loop is given by 

-g\TXaAT%b, j d'^yd^'z [Sf{x - y)^MT:^' G^,{y ~ z) 

X [5f(x-z)7,q(z)]^^^^ (4.59) 
where the free quark and gauge propagators are given in momentum space as 



The above contribution can be written as 
.2 r ^d^jD 



g^^ - (1 - A)— ^ 



(4.60) 



(4.61) 



where 



j^-j^ [Spip + fc)7p]cai Gf,^{k) [Spiq - fc)7i^]/3^i , (4-62) 
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whose divergent part is independent of the momenta p, q and is given by 



Taai,l3l3i (0, 0) 



1 1 



167r2 e 



XI ® 1 



with a^^i, = ^ [7/i,7i/]- The divergent part of the 1-loop eontribution becomes 



9 



,-[(To + ATi).g'^(a;)®g''(.T)]f (4.64) 



where (bold-faced symbols represent matrices in flavor and spinor space) 



fj,iy 



{i<E>i + Nmi}"'-Z • 



Here we use the notation 



By the following Fierz identities for spinor indices 

~ i ^A"^ ® c^^,. = Pr®Pr + Pl®Pl- 2{Pr®Pr + Pl®Pl) 
-\Y^ = Pr®Pr + Pl®Pl - 2(P^ «) Pr + Pl® Pl) 



(4.63) 



(4.65) 
(4.66) 

(4.67) 
(4.68) 

(4.69) 
(4.70) 



where Pr, Pl are the chiral projectors, Pr, = (1 + 75)72 and Pl = (1 — 75)72, Tq can 
be simplified as 

(4.71) 

where either ai,/3i G {1, 2} (right- handed) or ai,/3i G {3, 4} (left-handed) due to the 
chiral projections in eqs. (|4.69p and (|4.70p . In the following calculation of the 1-loop 
anomalous dimensions, eq. (|4.64p together with eqs. (|4.7ip and (|4.66p are the key 
equations. 

We now calculate the anomalous dimensions of general 3-quark operators at 1-loop. 
In terms of the renormalization factor defined as 



B. 



3 ■ ^Z^,[qt,]^Z^,Zy\q% Z3, = l+g^(Z^;^ + Z^;;j + ... , (4.72) 



.2/7(1) 



.(1) 



where Z^^ ( Z^^^^ ) is the A-independent (dependent) part at 1-loop, the divergent 

part of the insertion of the 3~quark operator = B^^^ at 1-loop is given by a linear 
combination of insertion of baryon operators as 
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(r(i)d'v)^ = -g' Zil^ + + '-Z^^'^ Bf: . (4.73) 



FF' 

, •^'7(1)1 

/ rr' 

The A-dependent contribution from Ti in (j4.66p is diagonal and given by 

.g^(rr^)F = 3A3f^^Bf, (4.74) 

so that the A-dependent part of vanishes: 

(1) __^N+1 3 (1) A 3(iV+l)(jV-3) 

327r2 TVe "2^^ ~647V7r2 e - (A^ - 3). (4.75) 

Therefore is A-independent, as expected from the gauge invariance. We leave N 
explicit in some formulae to keep track of the origin of the various terms, but we should 
always set iV = 3 at the end. 

The A-indepcndcnt part of r(i) from To in (|iTT|) leads to {N = 3): 

/p(l)divN/g'i _ (iV +1) g \„r^fgh r^i^fgh r^i^fgh or/s'i] (a 70\ 

2 



(Yii)<iiv^fgh ^ [N + 1) g 



(4.77) 



where a, /3,7 € {1,2} (right-handed), while 7 e {1 = 3,2 = 4} (left-handed). Note 
that the same results hold with hatted and unhatted indices exchanged. These results 
can be easily diagonalized as 



V^3g ){aal3} 


— fy^-'-h^s'' — 


e 


(4.78) 


(yC^)\f^gh 
\^3q )[al3]a 




-ni. 


(4.79) 


(yW\fgh 

y^3q ){al3}^ 


— ( y'-^h-f^ii _ 


4 


(4.80) 


(yW^fy^gh 
y^3q )[al3]j 


_ (^W\f7^gh _ 

^ 3(j '[ap]j 


-ni. 

e 


(4.81) 



where d is given by 

^"32^-9^- ^'-''^ 
The square bracket denotes antisymmctrization [a/3] = af3 — /3a, and curly bracket 
means {a/3} = a/3 + /3a, {aa/3} = aa/3 + a/3a + /3aa. The totally symmetric case 
corresponds to the decuplet representation (for Nf = 3) and contains the Nf = 2,1 = 
3/2 representation. The antisymmetric case corresponds to the octet representation 
(for Nf = 3) and contains the Nf ~ 2,1 = 1/2 representation. The anomalous 
dimension at 1-loop, obtained from 

7 = g'l^'^ + Oig') = /3Dig, ^)^^ = -^zi'Ja'e + 0(/) , (4.83) 

becomes 

(7^'^^'''' - (7^'^^'''' =-24d, (4.84) 

V J {aa 13} V J {aal3} 
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(7(1)) = =24d, (4.85) 



i/3]q V / [&l3]c 
Jgh / ,,.\fgh 



(^^^A ' = (j^^A ' = -Sd, (4.86) 
(t^'M =(7^') - =24rf. (4.87) 



We next consider the renormalization of arbitrary local gauge invariant 6-quark 
operator of (lowest) dimension 9, which can be written as a linear combination of 
operators 

Ocix) = B^lflix) = B^l{x)B^l{x) - Oa{x)Ob{x) , (4.88) 

with A ~ (ri,i^i) and B = (r2,F2). Note Oa{x) and/or Ob{x) may not be opera- 
tors with proton or nucleon quantum numbers and separately may not be diagonally 
renormalizable at one loop. The reason for considering the renormalization in more 
generality is that in principle there may be operators in this class which occur in the 
OPE of two nucleon operators at higher order in PT, but are relevant in the analysis 
because of their potentially large anomalous dimensions. 

According to the structure of the divergent part at 1-loop order, the operators in 
eq. (|4.88p mix only with operators Oc' ~ Oa'Ob' which preserve the set of flavors 
and Dirac indices in the chiral basis i.e. 

FiUF2= F{ U , Fi U Fa = F'^ U F'2 . 

Note however that such operators are not all linearly independent. Relations between 
them follow from a general identity satisfied by the totally antisymmetric cpsilon 
symbol which for N labels reads 

jYgai...ajVgbl...6jv _ ^ ^ ^ai ...aj_i !)fc a^ + i . . .ajv gf>l ■■■bfc- lOj 'Jfc+l ■ ■ ■&« ^/^ gg^ 

For our special case, = 3, this identity implies the following identities among the 
6-quark operators 

3i^S:S+E^EH)K=0' (4-90) 

where i-th index of abc and j-th index of def are interchanged in (abc,def)[i, j]. For 
example, (Fi,F2)[l,l] = a2/3i7i> "1/3272 or (Fi,F2)[2,l] aia27i, /3i/3272- Note that 
the interchange of indices occurs simultaneously for both Fi, F2 and i^i, F2 in the above 
formula. The plus sign in (j4.90p appears because the quark fields arc Grassmann. 

As an example of identities, let us consider the case that Fi, F2 = aa/3, a/3/3 {a ^ (3 
and Fi,F2^ ffg, ffg (/ 7^ g). The constraint gives 

nnffgjfg I (o r)\r>ffgjfg , nfffjgg , to ^ ^\uf3g,fff 

aaa.BBB ^ \ / aBB.aaB 



_ Ar>ffg,ffg I nfffjgg 1 Tjfgg,fff _ n (a Q^\ 



where minus signs in the first line come from the property that B^^'^^ — —B^^'^\ 
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An immediate consequence of the identity is that the divergent part of the A- 
dependent contributions, calculated from Ti in (|4.64[) . must vanish, after the sum- 
mation over the 9 different contributions from quark pairs on the different baryonic 
parts A, B is taken. The A-dependent part of the contribution of quark contractions 
on the same baryonic parts is compensated by the quark field renormalization. Thus 
the renormalization of the bare 6-quark operator is A-independent as expected from 
gauge invariance. 

We thus need only to calculate the contributions from To, which can be classified 
into the following 4 different combinations for a pair of two indices: 



where / g and a 7^ /3 e (1, 2) (Right) or e (3, 4) (Left). 

The computation can be made according to the following steps: 

i. ) Select the total flavor content e.g. 3/ + 3g or 4/ + 2g (/ ^5). These are the 
only cases for baryon operators with Nf = 2, but the approach is also applicable to 
more general cases {Nf > 2). 

ii. ) Given a flavor content classify all the possible sets of Dirac labels in the chiral 
basis e.g. 111223, 112234, ... It is obvious from the rules above that some have equiva- 
lent renormalization at 1-loop e.g. 111223 and 112223 with 1 o 2, and also those with 
hatted and unhatted indices exchanged e.g. 111223 and 133344. 

iii. ) For given flavor and Dirac sets generate all possible operators. Then generate 
all gauge identities between them and determine a maximally independent set {Si}. 

iv. ) Compute the divergent parts of the members of the independent basis: 



V.) Finally compute the eigenvalues and corresponding eigenvectors of 7"^ to de- 
termine the operators which renormalize diagonally at 1-loop. 

Some of the steps are quite tedious if carried out by hand. e.g. in the case 3/ -t- 3g 
and Dirac indices 112234 there are initially 68 operators in step iii.) with 38 indepen- 
dent gauge identities, and hence an independent basis of 30 operators. However all 
the steps above can be easily implemented in an algebraic computer program using 
MATHEMATICA or MAPLE. 

For Nf = 2 the quark f{g) has I3 = 1/2 ( —1/2). If an eigenvalue is non-degenerate 
the corresponding eigenvector belongs to a certain representation of the isospin group. 
If the eigenvalue is degenerate then linear combinations of them belong to definite 
representations. For the 3/ -|- 3g case they can have 7 = 0,1, 2, 3. Eigenvectors with 
/ = 0, 2 are odd under the interchange f g and those with 7 = 1,3 are even. The 
operators in the case 4/ + 2g have I3 ^ 1 and hence have 7 = 1,2, 3. The eigenvectors 
in this case can be obtained from those of the 3/ + 3g case by applying the isospin 
raising operator. 




(4.92) 
(4.93) 
(4.94) 




(4.96) 
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The complete list of eigenvalues and possible isospins arc given in Tables 14.1114.31 
The most important results are summarized as follows. 

1) For the 3/ + 3g (and 4/ + 2g) cases all eigenvalues 7^ < 48(i — i-f^, where 7Ar 
is the 1-loop anomalous dimension of the nucleon (3-quark) operator. 

2) It is easy to construct eigenvectors with eigenvalue 2^n e.g. operators of the 
form B^fa tB^^( - since there is no contribution from diagrams where the gluon line 

joins quarks in the different baryonic parts. 

3) Operators with higher isospin generally have smaller eigenvalues. 

We now consider the rcnormalization of 6-quark operators which appear at the 
tree level of the OPE more explicitly. Since 



C75 



/O -1 \ 
10 
-1 

\0 1 / 



in the chiral representation, the nucleon operator is written as 



Bl^B 



ffg 

■-+a.[2.1] 



B 



ffg 

ct+aj2,i] 



(4.97) 



(4.98) 



where a = 1, 2, d = a + 2, and f ^ g. This has anomalous dimension 77V = 24d. 

Two independent 6-quark operators occurring in the OPE at tree level can be 
decomposed as follows. The spin-singlet {S = 0) and isospin-triplct (/ = 1) operator 
becomes 



B 



ffg 



B 



ffg 



a+a,[l3,a] + [P,&] P+P,[l3,a] + [P,&] 



B 



II 



B 



III 



B 



IV 



B'i 



B 



VI 



(4.99) 



where {a ^ /3) 



Bf 


_ j^f fg 


Triffg 


VB{f? 


nffg 

WAV 










(4.100) 


Bll 


_ j^f fg 


j^ffg 


f B^fl. 


nffg 
"WM ~ 


f B^^? 


nffg 


a[l3,a 


r>ffg 

mA 


,(4.101) 


Bill 


_ j^f fg 


r>f fg 


I nffg 

^ &lp.a 


T^ffg 

"mM ' 


f B^f? 


nffg 

mA] 


f B^f? 


f>ffg 

mA 


,(4.102) 


Biv 


_ ]^f fg 


r>ffg 


I nffg 

^ alP.a 


T^ffg 

WM ' 










(4.103) 


rOi 

By 


_ jjf fg 

Ct[/3,Q 


f)ffg 


I nffg 

^ a[P,a 


T^ffg 

P[PA\ ' 










(4.104) 


rOi 

Byi 


_ jjffg 

Q[/3,a 


nffg 


f i?^f|. 


T^ffg 
"WM ■ 










(4.105) 



In the above some contributions arc obtained from interchanges under (1,2) -H- (3,4) 
or (1,3)0(2,4). 

Similarly the spin-triplet {S — 1) and isospin-singlet (/ = 0) operator is decom- 
posed as 



B 



ffg 



B 



ggf 



a+&,[l3,a] + [P,a] a+a,[l3 ,a] + [l3,a] 



= B 



B 



II 



B 



III 



B 



IV 



Bi 



B 



VI ' 
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(4.106) 

where 







Tjggf 

a[P,a] ' 


f B^f|. 

a[p,a 


nggf 










(4.107) 


nlO 
BlI 


= -B^/l 

a[p.a 


Q:[/3,d] 


f s^f|_ 

a[p,Q 


q[P ,aj 




nggf 


^ a[/3,a 


Tjggf 


,(4.108) 


nlO 
^III 


= B^m 

a[p,a 


a[p.a\ 


1 n//9 

^ a[/3,a 


B''J. 

q[P ,aj 


f B^^f 

a[/3,Q 


j^ggf 


f B^^? 


j^ggf 

&[$.& 


,(4.109) 


rIO 


— B^^? 

a[/i?,a 


B'H - 


, n//9 












(4.110) 


nlO 

By 


_ fg 


B'U 1- 


, n//9 


j^ggf 










(4.111) 


nlO 




Q[/3,a] 


f s^/l. 

a[p,a 


T^ggf 

a[l3.a] ■ 










(4.112) 



It is important to note here that operators By^j for both cases {SI = 01 and 10) 
have the maximal anomalous dimension at 1-loop, since as noted in point 2) above, 
no 1-loop correction from Tq joining quarks from the two baryonic components exists 
for type of operators. Therefore there always exist some operators with 

f3A — which dominate in the OPE at short distance. 

The 1-loop corrections F^^^ to 6-quark operators B^^ are summarized as: 

(rn^'^ = -i2V, (r?})<^^ = i2^B?l, (r?},)'^' = o, 

(r?^)^'' = o, (r°i)('' = 6^< + 6^s°i,, (roi,)W = 24^i?ov, (4.113) 

for 5/ = 01. The last two results can be written as 

{n\f^ =6-^B'y\, (r«,V)''^=24^i?0V,, (4.114) 

where 

B"y\^B"y'-lB"y'j, B^^, ^ B^^ . (4.115) 

Similarly for SI = 10 

[Vfyf^ ^%'iBf^, (rio,)^'^=6-i?iP,, (ri?,,)*'^=24-<,,, (4.116) 
e e e 

where 

Bi", = < - , B},\ = B]^, . (4.117) 

Denoting the eigenvalues of the anomalous dimension matrix by 7c, the values of 
7"^^ defined by 
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7c - 2jN = 2d7^ , (4.118) 

are given in table (iV = 3), which shows, in both cases, that the largest value is 
zero while others are all negative. The case 2 in sect. 14.11 is realized: f3c = and 

/3c' =/3o"' = -^^7^ for 5 = 0,7=1, (4.119) 
/3c'=/5o'" = -^^-^ for 5= 1,7 = 0. (4.120) 

4.4 Short distance behavior of the potentials and the repulsive 
core 

As discussed before, the NN potential at the leading order of the derivative expansion 
is given by 

V\r) = V„'ir) + l//(r)ai • + Vi{r)Si2 + 0(V). (4.121) 

From the result in the previous subsection, the OPE of NN at tree level can be written 

as 

Bl{x + y/2)Bl{x ~ y/2) ~ cv/73(«,_„^(x) + cn{- \ogrf"' B{l^p{x) + • {4.122) 

where cyi and c// are some constants, and • • • represents other contributions, which 
are less singular than the first two at short distance. The spinor and flavor indices 
a, P and /, g are explicitly written here. The anomalous dimensions /3q ^ are given in 

Kim and mm . 

In the case that S = Q and 7 = 1 (a 7^ /3 and / = g in eq. (|4.122p ). the leading 
contributions couple only to the J = L = state given by 

\E) = \L, = 0, 5, = 0, h = 1)l=o,s=oj=i - |0, 0, l)o,o,i • (4.123) 

The relevant matrix elements are 

cv/(0|73^1„^|0, 0,1)0,0,0 = AljYo^[a(3]{fg)^, (4.124) 
c//(0|B/f,„/3|0,0, 1)0,0,0 = A?,rooM{/<?}i, (4.125) 

where and are non-perturbative constants, Yll^ is a spherical harmonic 

function, [a/3] = {5ai5i32 — '5;3i<5q2)/a/2 represents the [S, Sz) = (0,0) component, and 
{/.9}i = 5fi5gi corresponds to isospin {I,Iz) ~ (Ijl)- The wave function at short 
distance is dominated by 

'PE^'iy) = {0\Bi{x + y/2)B^x~y/2)\0,0,l)o,oA 
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:^(A^ + <(-logr)ft';^)0(i5o,O)" + --., 
where 0(^50, Jz = 0)^^^"^^ — loo[a/?]{/.9}i- This wave function leads to 



II j.n 



rriN 



(4.126) 



(4.127) 



where m = mpj /2 is the reduced mass of the two nucleon system. Since 6*12 is zero on 
(^(i^cO)", we have 



where v (r) — Vo(r) — 314^ (r). Therefore the potential is obtained as 



where 



-/3o^^(-logr)^o'-i 



(4.128) 



(4.129) 



(4.130) 



The potential diverges as i^°^(r) in the r — J- limit, which is a little weaker than r~^. 

In the case of the spin-triplet and isospin-singlct state (S' = 1 and / = 0), the 
leading contributions in ea. (|4.122p couple only to the J ~ 1 state, which is given by 



1^^) = p5i,J, = l)+x|3Z3i,J, = 1) 



where 



l^^i, J, = 1) = \L, =0,Sz = 1)l=o,s=i , 
1 



|0,1) - V3|l,0) +\/6|2,-l) 



L=2,S=1 



(4.131) 

(4.132) 
(4.133) 



X is the mixing coefficient, which is determined by QCD dynamics, and '^^'^^Lj specifies 
quantum numbers of the state. We here drop indices I and Iz unless necessary. 
Relevant matrix elements are given by 



B 



fa 

i.af3 



B, 



fg 



Di,l 



for i ~ II and VI, where Bf are non-perturbative constants, and 

0('5i) =yoo(e,0){a/3}i[/g], 

with {a/3}i = 6aid/3i. 



0, 



(4.134) 



(4.135) 
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Using the above results, the BS wave function becomes 



which leads to 



i(y)^{B°., + (-logr)^"i?°,}^e50, 
-/310 (-logr)^""-! 



(4.136) 



2m 



rriN 



^ir^y) c ^ ^ Sy,^ (-^Si) . (4.137) 

(4.138) 

Y2o{aP}i-VSY2i{aP}o + V6Y22{al3}^i [fg], (4.139) 



On the other hand, 
where T/i°(r) = V^{r) + V°{r) and 



with {a/3}_i = Sa2Sp 2 and {a/3} p = (^qi££2 + 1(502 )/\/2. 
By comparing ea. (|4.137p with ea. (|4.138p . we obtain 



Ki°(r)^Fi°(r)^ 
^vi 



VTir)^0. (4.140) 

This shows that the central potential V^^{r) diverges as F^'^{r) in the r — >■ limit, 
which is a little weaker than 1/r^, while the tensor potential Vrir) becomes zero in 
this limit at the tree level in the OPE. 

While the OPE predicts the functional form of the central potentials at short 
distance fro both 5* = and 5 = 1 channels, it can not determine whether it is 
repulsive or attractive. The NN operators are decomposed as 



BiBf, = [Bi + Bii + Bin + Biv + By + B 



vi_ 



(4.141) 



where 





~ [BbrBrr - 


^ BllBll]^^!} , 






(4.142) 




= [BrrBrl - 


h BrlBrr + BllBlr - 


h BlrBll 


fg 

a/3 ' 


(4.143) 


iBiii)% 


= [BrrBlr - 


V BlrBrr + BllBrl - 


\- BrlBll 


fg 

ap 1 


(4.144) 


{Biv)% 


~ [BrlBrl H 


- BlrBlrYo/^ , 






(4.145) 


{Bv)% 


= [BrlBlr H 


-BlrBrl\% , 






(4.146) 


(BviYj, 


= [BrrBll H 


- BllBrr^J'^ . 






(4.147) 



Here 3-quark operators in terms of left- and right- handed component are defined by 



^AX,a ~ 



{PaUB^p , = Bif^{C^,Px)p^{iT2)sh , 



fgh . 



for A, X ^ R 01 L. We need to know 

{0m)i$\2N,E) 

for i = II, VI. 



(4.148) 
(4.149) 
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For f ^ g, Lorentz covariance leads to 

(0|S^B^|2N,£) = C^Y{s)PAu{p,ai)PBu{-p,<j2), (4.150) 

A,B=R,L 

where s ~ E"^ ^ 4(p^ + m^) with the total energy E in the center of mass frame, 
(Ti {i = 1,2) is the spin of the i-th nucleon, and C^y is an unknown function of s. 
Note that spinor indices are suppressed here. Invariance of QCD under the parity 
transformation PBxP~^ ~ liBx where R ^ L and L = R gives 

(|4:T501) = {0\PBf,B^p-'P\2N,E) ^Y.^xfPAlM-P,<^i)PBlMP,'^2) 

A.B 

= IZ^if^A^i(p,^i)PB^i(-p,^2), (4.151) 



where 74u(— p, ui) = u(p, ai) is used. The above relation implies = C^y - There- 
fore the matrix elements are evaluated as 



(0|(S//)^f '^|2N, E) = C'^t+LRiiPR ® Pb +Pl® ^l)"(P, ^l)«(-P, ^2)}a/3T^a 



(4.152) 



and 



{'^\{Bvi)if'^m.E) = C'^l' {{Pr ®Pl + Pl® PR)u{p,ai)u{-p,<j2)}c.PTPo. 

(4.153) 

The spinor in Dirac representation for 7 matrices ptzykson and Zuber, 1980] ) is 
given by 



it(±p, 



1 



y/Eis! + rriN 



/ En + rriN \ 

TPz 

V / 



,m(±P, -) = 



1 



V -Eat + mN 



( \ 

En + rriN 



V J 

(4.154) 



for p = (0,0,p3 > 0), where En .Jp^ +m%. For / = 1 ( fg + gf) and 5* = 
(di = + and cr2 = — ) the above explicit form for the spinors gives 



{{Pr ® Pr + Pl® PlMp, +)u{-p, -)}i2--21 = En , 
{{Pr (S)Pl + Pl® Pr)u{p, +)u{-p, -)}i2-2i = ruN 



(4.155) 
(4.156) 



while, for / = ( fg — gf) and 5=1 (cti = + and 0-2 = + ) 



{{Pr ®Pr + Pl® PlHp, +)u{-p, +)}n = mN , 



(4.157) 
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{{Pr. ®Pl + Pl® PrHp, +)u{-p, +)}ii = En ■ 



(4.158) 



Finally the ratio of the matrix elements becomes 




(4.159) 



for fg + gf and (cri, cr2) = (+, -) ( 



^So ), and 



(0|(i?,,){f-^^|2N,i^) _ niN C^RL+LRis) 



(4.160) 



for fg - gf and (cti, 0-2) = (+, +) ( ^Si ), where s = 4E%. 

Unfortunately, the sign of the ratio for these matrix elements can not be deter- 
mined. As a very crude estimation, the non-relativistic expansion for constituent 
quarks whose mass mg is given by ttiq = totv/S is considered. In the large mg 
limit, 74(70 = Qo and 74U0 — "Oi where a subscript for q and u means the 0-th 
order in the non-relativistic expansion. In this limit, it is easy to show C^y = ^ 
for all X,Y,A,B, so that Cff^.^^ = 2C and Cff = C. Furthermore the first order 
correction to Cxy = C* vanishes in the expansion. Therefore the leading order of the 
non-relativistic expansion gives 




(4.161) 



for (cti,CT2) = (+,-) ( 



^5*0 ), and 




(4.162) 



for ((Ti,CT2) = (+, +) ( '^5'i ). For both cases, the ration has positive sign, which gives 
repulsion at short distance, the repulsive core. 
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Table 4.1 Eigenvalues 7j of the anomalous dimension 
matrix 7 and isospins of the corresponding eigenvec- 
tors for the case 3f3g. 



Dirac indices 


7,/ (2d) 


/ 


mill 


-24 





111112 


-24 


0,1 


111122 


-4 







-24 


0,1 




-40 


2 


111222 


-4 







-12 


1 




-24 


0,1 




-40 


2 




-72 


3 


111113 


-16 


0,1 


111123 


-6 


0,1 




-16 


0,1 




-24 


1,2 


111223 





0,1 




-6 


0,1 




-16 


0,1 




-18 


1,2 




-24 


1,2 




-48 


2,3 


111133 


-4 







-16 


0,1,2 


111134 





1 




-4 







-12 


1 




-16 


0,1,2 


111233 


4 


1 




-4 







-8 


0,1,1 




-16 


0,1,2 




-32 


1,2,3 
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Table 4.2 As in TabljO (continued!. 



Dirac indices 


7j7(2d) 


/ 


111234 


20 







8 


1 




4 


1 







1 




-4 







_s 
o 


119 

U, i, i, z 




— 12 


1 




-16 


0, 0,0, 1, 1, 2, 2, 2 




-32 


1,2,3 


112233 


8 







4 


1 




-4 


0,0,1,2 




-8 


0,1,1,2 




-16 


0, 1,2 




-28 


2 




-30 


1,2,3 


112234 


20 







1 9 

IZj 


1 




8 


0,1 




A 


1 
± 







1,1 




-4 


0,0,1,2 




-8 


0,1,1,2 




-12 


1 




-16 


0,0,1,1,2,2,2 




-28 


2 




-32 


1,2,3 




-36 


1,2,3 


111333 


-6 


0,1 




-24 


0,1,2,3 


111334 





0,1,1,2 




-6 


0,1 




-18 


1,2 




-24 


0,1,2,3 
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Table 4.3 As in TablegT] (continued). 



Dirac indices 


7,7 (2d) 


/ 


112334 


24 


0,1 




6 


0,1 







0,0,1,1,1,1,2,2 




-6 


0,1 




-12 


1,1,2,2 




-18 


1,1,2,2 




-24 


0,1,2,3 




-30 


0,1,2,3 



Table 4.4 The value of 7^-^ (defined in (4.72)) for 
each eigen operator in the SI = 01 and SI — 10 states. 





/ 


// 


/// 


IV 


V 


vr 




-36 


-12 


-24 


-24 


-18 





^10 


-28 


-4 


-24 


-16 


-18 






5 

Concluding Remarks 



In this lecture, we consider two different but closely related topics in lattice QCD 
approaches to hadron interactions. 

In the first part, we explain the Liischer's formula, which relates the scattering 
phase shift to the two particle energy in the finite volume. Although this formula 
is well-established and has appeared in the well- written reference ( jLiischer , 1 99 1 [ ) , a 
comprehensive (but less rigorous) derivation for the formula has been attempted in this 
lecture for the tttt system as an example with the emphasis on the BS wave function. 
It is important to stress that the BS wave function at large separation behaves as 
the free scattering wave with the phase shift which is determined by the unitarity of 
the 5'-matrix in QCD. The Liischer's formula can be obtained from this asymptotic 
behavior of the BS wave function. 

In the second part, on the other hand, we consider the BS wave function at non- 
asymptotic region where the interaction between two particles exists, in order to define 
the potential in quantum field theories. We apply this method for the two nucleons 
to calculate the NN potential in QCD. The first result from lattice QCD has a good 
shape, which reproduces both the repulsive core at short distance and the attractive 
well at medium and long distances. 

In the last part, the origin of the repulsive core is theoretically investigated by the 
OPE and the renormalization group. The analysis predicts the r dependence of the 
potentials at short distance, though it can not tell the sign of the potential, positive 
(repulsive) or negative (attractive). A crude estimate by the non-relativistic quark 
gives the positive sign, the repulsive core. 

The method to investigate potentials through the BS wave function in lattice QCD 
is a quite new approach, which may be applied to the following directions. 

(1) In order to extract the realstic NN potentials from lattice QCD, it is necessary 
to carry out full QCD simulations near the physical u, d quark masses. Studies along 
this line using 2-|-l-flavor QCD configuration generated by PACS-CS Collaboration 
( |Aoki et ai, 2009a| is currently under way( |Ishii et ai, 2008[ ). 

(2) The hyperon-nucleon {YN) and hyperon-hyperon (YY) potentials are essential 
to understand properties of hyper nuclei and the hyperonic matter inside the neutron 
star. While experimental scattering data are very limited due to the short life-time 
of hyperons, the NN, YN and YY interactions on the lattice can be investigated 
in the same manner just by changing only the quark flavors. The 5iV potential in 
quenched QCD( Nemura et ai, 2009| and the AA'' potentila in both quenched and full 
QCD(Nemura et ai, 20081 are examined as a first step toward the systematic un- 
derstanding of baryonic potentials. To this connection, the OPE analysis should be 
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extended to the Nf = 3 case, in order to reveal the nature of the repulsive core in 
baryon-baryon potentials. Since quark mass terms can be neglected in the OPE at 
short distance, the analysis can be done in the exact SU(3) symmetric limit. 

(3) The three-nucleon force is thought to play important roles in nuclear structures 
and in the equation of state in high density matter. Since the experimental informa- 
tion is very limited, the extension of the method to the three nuclcons may lead to the 
first principle extractions of the three-nucleon potentials in QCD. It is also interesting 
to investigate the existence or the absence of the repulsive core in the three-nucleon 
potentials. The calculation of anomalous dimensions of 9-quark operators will be re- 
quired at 2-loop level. 

(4) More precise evaluations including numerical simulations of matrix elements 
(0|Oa'|^) will be needed to understand the nature of the core in the potential through 
the OPE analysis. 
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